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I.  INTRODUCTION  AND  SUMMARY 


For  a  theoretical  description  of  power  plants  that  derive  their  energy  from  the  combustion  of  a 
fuel,  it  is  necessary  to  solve,  to  an  appropriate  degree  of  precision,  a  hydrodynamic  problem  that 
requires  for  its  solution  a  knowledge  of  the  thermodynamic  properties  of  the  working  fluid  com¬ 
posed  of  the  products  of  the  combustion  reaction.  It  is  usually  a  good  approximation  to  assume 
that  the  properties  of  the  combustion  gas  are  determined  by  the  conditions  of  thermal  equilibrium, 
and  it  is  thus  possible  to  employ  the  methods  of  classical  thermodynamics  for  their  computation. 
The  thermodynamic  properties  of  fuel  gases  also  are  of  considerable  importance,  as  they  form  the 
basis  for  the  design  of  appropriate  means  for  their  effective  utilization.  The  specification  of  the 
operating  conditions  to  produce  a  gas  for  use  as  an  intermediate  in  a  chemical  process,  such  as 
the  synthesis  of  liquid  fuels,  may  be  based  upon  a  study  of  the  variation  of  the  composition  of 
the  synthesis  gas  with  changes  in  the  various  process  variables. 

Although  there  exist  a  large  number  of  important  scientific  and  technical  applications  of  the 
data  that  can  be  obtained  from  the  systematic  determination  of  the  thermodynamic  properties  of 
combustion  gases,  such  application  has  been  greatly  handicapped  by  the  extremely  tedious  and 
time-consuming  computational  methods  required.  The  development  of  large-scale  automatic  com¬ 
putational  equipment  makes  feasible  the  initiation  of  a  systematic  program  for  the  determination 
of  the  thermodynamic  properties  of  combustion  gases  and  for  the  application  of  such  data  to  spe¬ 
cific  problems  of  scientific  and  technical  importance. 

As  a  part  of  its  program  of  basic  research  in  flame  and  combustion  phenomena,  the  Explo¬ 
sives  and  Physical  Sciences  Division  of  the  Bureau  of  Mines  is  engaged  in  such  a  program.  The 
present  communication  is  designed  to  introduce  a  series  of  reports  presenting  the  results  of  this 
program. 

There  has  long  existed  a  need  for  systematic  and  economical  methods  for  the  calculation  of 
the  thermodynamic  properties  of  systems  of  many  constituents,  and  this  need  is  emphasized  by 
the  application  of  automatic  computational  equipment.  In  the  present  report,  we  wish  to  assemble 
the  methods  that  have  been  employed  in  this  laboratory  and  that  have  been  shown  to  be  particu¬ 
larly  well-adapted  to  routine  application. 

Consideration  of  the  thermodynamics  of  combustion  gases  is  introduced  by  a  discussion  of 
the  stoichiometry  of  multicomponent  systems  and  of  applications  of  the  phase  rule  to  such  sys¬ 
tems.  This  formal  treatment  serves  to  introduce  the  notation  that  will  be  employed  consistently 
throughout  the  remainder  of  the  report,  and  it  establishes  the  conditions  for  equilibrium  in  multi- 
component  systems  in  a  mathematical  form  that  is  particularly  appropriate  for  computational 
procedures. 

Two  methods  are  described  for  the  calculation  of  the  equilibrium  composition  of  systems  of 
many  constituents.  These  methods  are  employed  in  the  formulation  of  several  specific  compu¬ 
tational  programs  for  the  calculation  of  the  equilibrium  composition  of  systems  composed  of  com¬ 
pounds  of  carbon,  hydrogen,  oxygen,  and  nitrogen. 


Throughout  this  report,  we  assume  that  the  ideal  gas  law  provides  an  adequate  equation  of 
state  for  the  gaseous  phase.  However,  it  is  shown  that  the  methods  for  the  calculation  of  equi¬ 
librium  composition  are  easily  generalized  to  take  account  of  gas  imperfections,  and  that  this 
generalization  does  not  affect  in  any  important  respect  the  detailed  computational  programs  al¬ 
ready  formulated.  The  necessary  general  relations  are  summarized.  They  can  be  applied  to  spe¬ 
cific  forms  of  the  equation  of  state  without  difficulty. 

The  calculation  of  the  thermodynamic  properties  of  the  equilibrium  mixture  is  straightforward 
after  the  composition  of  the  mixture  is  known.  We  summarize  briefly  the  well-known  thermody¬ 
namic  relations  that  are  employed  for  this  purpose. 

We  conclude  with  a  description  of  a  massive  generalized  table  of  the  equilibrium  composition 
and  thermodynamic  properties  at  equilibrium  of  the  four- component  system  containing  compounds 
of  carbon,  hydrogen,  oxygen,  and  nitrogen.  The  table  will  list  the  composition  and  thermody¬ 
namic  properties  over  an  extended  range  of  temperature  and  pressure,  and  it  will  include  all  pos¬ 
sible  combinations  of  the  four  elements,  the  spacing  of  calculated  points  being  small  enough  to 
permit  accurate  interpolation.  The  table  will  be  the  result  of  about  half  a  million  individual  com¬ 
putations,  and  its  construction  utilizes  the  computational  methods  described  in  this  report. 
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III.  STOICHIOMETRY  AND  APPLICATION  OF  THE  PHASE  RULE 


The  composition  and  the  thermodynamic  properties  of  a  mixture  at  equilibrium  are  independent 
of  the  path  by  which  equilibrium  is  attained.  From  the  point  of  view  of  thermodynamics,  the  mo¬ 
lecular  form  in  which  the  different  elements  are  initially  introduced  to  the  system  is  therefore  ir¬ 
relevant,  and  it  suffices  to  specify  the  gross  composition  of  the  system  in  terms  of  the  elements 
contained  by  it.  As  the  composition  of  the  system  at  equilibrium  is  most  appropriately  expressed 
in  terms  of  intensive  quantities  (mole  fractions),  the  results  are  independent  of  the  total  size  of 
the  system,  which  may  be  taken  to  be  any  convenient  value.  Therefore,  if  the  number  of  different 
elements  in  the  system  is  m,  it  is  sufficient  for  a  unique  description  of  the  system  to  specify  the 
number  of  gram  atoms  of  m-1  elements  relative  to  an  arbitrary  amount  of  the  m-th  element,  i.e.,  to 
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specify  values  of  the  quantities  (x(k)/x(m))  k  =  1,  2,  ...  m-1,  where 


* 


Qk  =  (X(k)/X(m>)  Qm)  (3.1) 


and  and  Qm  are  the  number  of  gram  atoms  of  the  k-th  and  m-th  elements,  respectively,  avail¬ 
able  to  the  system.  Alternatively,  one  may  specify  m-1  different  linearly  independent  combina¬ 
tions  of  these  ratios  Rkm.  We  note  that  it  is  not  always  necessary  to  specify  as  many  as  m-1  dif¬ 
ferent  quantities. 

Description  of  Carbon,  Hydrogen,  Oxygen,  Nitrogen  Systems 


A  large  number  of  different  sets  of  parameters  can  be  employed  for  the  description  of  the 
gross  composition  of  any  given  system.  In  this  section,  we  introduce  certain  parameters  that  are 
particularly  convenient  for  the  description  of  systems  composed  of  compounds  of  carbon,  hydro¬ 
gen,  oxygen,  and  nitrogen,  and  we  deduce  several  relations  that  will  be  employed  later.  Alterna¬ 
tive  descriptions  may  be  transposed  into  the  language  of  this  section  by  obvious  algebraic 
procedures. 

For  a  discussion  of  the  thermodynamics  of  combustion  gases,  it  is  convenient  for  the  descrip¬ 
tion  of  the  gross  composition  of  the  system  to  employ  the  parameters, 

A  =  Qn/(4Qc  +  Qh  +  Qn)> 

B  =  QH/(4QC  +  Qh>,  (3-2) 

C  =  2Q0/(4QC  +  QH  +  2Qo), 

where  Qq  Q^,  Qq,  Qjj  are  the  number  of  the  gram  atoms  of  carbon,  hydrogen,  oxygen,  and  nitro¬ 
gen,  respectively,  available  to  the  system.  Quantity  C  is  a  measure  of  the  amount  of  oxygen 
actually  available  to  the  system  relative  to  that  required  for  conversion  of  all  of  the  carbon  to 
carbon  dioxide  and  all  of  the  hydrogen  to  steam.  At  this  stoichiometric  point  C=  1/2.  Rich  fuel- 
oxidant  systems  are  characterized  by  C  <l/2  and  lean  systems  by  C  >1/2. 

It  will  be  shown  that  for  the  characterization  of  rich  fuel -oxidant  systems  consisting  of  a  gas 
mixture  in  heterogeneous  equilibrium  with  solid  carbon,  it  is  sometimes  convenient  to  describe 
the  relative  hydrogen,  oxygen,  and  nitrogen  content  by  means  of  parameters  that  are  independent 
of  the  amount  of  carbon.  For  this  purpose,  an  appropriate  set  of  parameters  consists  of 

b  =Qh/(4Qc  +  Qh) 

D  =  Qn/(  Qh  +  2Qo  +  On)  (3-3) 

E  =2Qo/(Qh  +  2Q0)- 

We  note  the  sets  (3.2)  and  (3.3)  are  not  independent,  but  that 
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D  = 


_ A(l-C) _ 

( 1-A)  [B(l-C)  +  C]  +  A(l-C) 


E  = 


C 

B(l-C)  +  C. 


-  3  - 


(3.4) 


Equations  (3.2)  or  (3-3)  determine  the  amounts  available  to  the  system  of  three  of  the  elements 
relative  to  the  amount  available  of  the  fourth.  If  oxygen  is  selected  as  the  reference  element, 
then 


(C/O) 


QC  ( 1-B)  Cl-C) 

Qo=  2C 


(1-B)  (1-E) 
2BE 


Qh  2B(  1-C)  2(  1-E) 

(H/°)  =  qq  =  q  =  £ 


(3.5) 


(N/O) 


Qn  2A(1-C)  2D 

Q0  =  0.1- A)  =  E(i-D) 


For  the  characterization  of  fuel-air  systems,  it  is  convenient  to  employ  the  parameters  A,  B, 
and  C  as  defined  by  equations  (3-2)  to  describe  the  relative  concentrations  of  carbon,  hydrogen, 
oxygen,  and  nitrogen  of  the  fuel  alone  without  regard  to  the  constituents  of  air.  These  fuel  pa¬ 
rameters  may  be  distinguished  from  their  more  general  counterparts  by  a  superscript  f.  Then,  if 
the  fuel-air  ratio  is  determined  by  the  general  parameter  C,  taking  account  of  the  total  available 
oxygen,  it  is  easy  to  show  that 


b(l-Bf)(l-C)(l-Cf)  +  c  (C-Cf) 

(C/O)  =  - - - 7 - 7 - > 

2bC(l-Cf)  +  2c  (CCf) 

4bBf(  1-C)  ( 1-Cf) 

(H/O)  =  - 7 - 7~ 

2bC(l-Cf)  +  2c(CCf) 

4bAf(l-C)  (1-Cf)/(1-Af)  +  2a(C-Cf) 

(N/O)  =  - - —7 - 7 - 

2bC(l-Cf)  +  2c(C-Cf) 

d  (C-Cf) 

(A/O)  =  - - 

2bC(  1-Cf)  +  2c  (C-Cf)’ 


(3.6) 


where  A^,  B^,  C^  are  measures  of  the  fuel  composition,  C  is  the  measure  of  the  total  amount  of 
available  oxygen  in  fuel  and  air,  and  where  a  =0.780881,  b  =0.209495,  c  =0.000300,  and  d  = 
0.009324.  The  composition  by  volume  of  dry  air  has  been  taken  to  be 

Oxygen,  20.9495%  Carbon  dioxide,  0.0300% 

Nitrogen,  78.0881%  Rare  gases  (Argon),  0.9324% 


If  the  fuel  is  a  hydrocarbon,  A^  -C^  -O,  and  equations  ( 3.6)  reduce  to 


(C/O) 


b(l-Bf)  (1-C)  +  cC 
2C  (b  +  c) 
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(H/O) 


4bBf(  1-C) 
2C(b  +  c)’ 


(N/O) 


a 


b  +  c 


,  (A/O) 


d 

2  (b  +  c) 


(3.7) 


The  Components 


In  a  system  containing  many  constituents,  it  is  possible  to  select  certain  constituents  that 
are  sufficient  for  the  complete  thermodynamic  description  of  the  composition.  By  this  it  is  meant 
that  if  the  system  is  conceived  to  consist  of  the  selected  constituents  only,  its  gross  composi¬ 
tion  (in  terms  of  die  amounts  of  each  chemical  element  present)  is  completely  defined.  The  con¬ 
stituents  thus  sufficient  to  describe  the  composition  are  called  the  components  of  the  sy  item. 
An  analytical  criterion  has  been  published-3/  for  the  choice  of  the  components.  In  terms  of  this 
criterion,  the  conditions  for  equilibrium  may  be  written  in  a  form  that  has  a  high  degree  of  sym¬ 
metry  and  is  particularly  well-adapted  for  formulating  a  computational  method  for  the  calculation 
of  the  equilibrium  composition. 

The  number  of  constituents  of  any  system  depends  upon  the  accuracy  with  which  it  is  desired 
to  describe  its  composition.  The  constituents  to  be  considered  must  be  chosen  a  priori,  and  this 
choice  usually  will  imply  the  neglect  of  certain  equilibria  that  may  be  expected  to  exert  a  negli¬ 
gible  effect  on  the  composition  of  the  system  at  equilibrium. 

Consider  a  closed  system  containing  s  different  substances,  which  are  assumed  to  be  in 
chemical  equilibrium.  The  molecular  formula  of  the  i-th  substance  may  be  represented  by 


Y(i)  =  yd)  .  .  •  X(k>  .  .  .  X<m).  , 
YW  XV'..  (Xik  «*"«’ 


ail 


aim 


(3.8) 


i  =  1,  2,  .  .  .  s,  where  X^k^  is  the  symbol  of  the  k-th  element,  a^  is  the  subscript  (which  may  be 
zero)  to  this  symbol  in  the  formula  to  the  i-th  substance,  and  m  is  the  total  number  of  elements 
represented  in  the  system.  For  every  i,  the  array  of  subscripts  a^,  k  =  1,  2,  .  .  .  ,  m,  may  be 
said  to  define  a  vector 


Yi  (®il>  •  •  ■  >  ®ik>  •  *  *  >  °"im^  (3-9) 

which  may  be  called  the  formula  vector  of  substance  i.  If  the  rank  of  the  matrix  of  the  vector 

elements  is  c,  it  follows  from  a  well-known  theorem  of  algebra  that  there  are  c  linearly  inde¬ 
pendent  vectors,  and  if  c  <s  there  are  (s-c)  linearly  dependent  vectors,  which  may  be  expressed 

as  linear  combinations  of  the  independent  vectors.  It  may  be  assumed  that  the  independent  vec¬ 
tors  are  designated  by  the  values  1,  2,  ...  ,  c  of  their  index.  Then  the  dependent  vectors  may 
be  expressed  as  linear  combinations  of  the  form 


3/  Brinkley,  S.  R.,  Jr.,  Note  on  the  Conditions  of  Equilibrium  for  Systems  of  Many  Constituents: 
Jour.  Chem.  Phys.,  vol.  14,  1946,  pp.  563-64. 
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c  (3.10) 

2l  vij  y\  =yi> 

j=i 


i  c  +  1,  c  +  2,  .  .  .  ,  s.  To  equation  (3.10)  there  correspond  (s-c)  conceivable  chemical 
reactions 


^  V j j  yO)  =Y(i), 

j=l 


(3.11) 


resulting  in  the  formation,  from  the  c  substances  with  linearly  independent  formula  vectors,  of  the 
(s-c)  substances  with  linearly  dependent  formula  vectors.  It  follows  that  the  specification  of  c 
substances  such  that  their  formula  vectors  are  linearly  independent  is  sufficient  for  a  description 
of  the  composition  of  the  system.  Therefore,  the  number  of  components  of  the  system  equals  the 
rank  c  of  the  matrix  of  the  subscripts  to  the  symbols  of  the  elements  in  the  formulas  of  the  sub¬ 
stances  comprising  the  system.  The  number  of  components  is  usually,  but  not  necessarily,  equal 
to  the  number  of  different  elements  in  the  system.  It  may  be  noted  that  the  choice  of  independent 
vectors  is  not,  in  general,  unique  and  that,  in  consequence,  the  choice  of  c  substances  as  com¬ 
ponents  and  the  expression  of  the  remaining  (s-c)  substances  as  products  of  reactions  involving 
only  the  chosen  components  is  usually  not  unique.  However,  for  the  heterogeneous  case,  con¬ 
sisting  of  a  gas  mixture  in  equilibrium  with  a  solid  phase  composed  of  a  single  constituent,  these 
criteria  can  be  satisfied  only  by  the  selection  of  that  constituent  as  a  component. 

This  discussion  has  demonstrated  the  possibility  of  a  choice  of  components,  which  make  it 
possible  to  express  each  of  the  dependent  constituents  as  products  of  reactions  involving  com¬ 
ponents  only.  Our  computational  procedure  for  the  determination  of  the  equilibrium  composition 
is  based  upon  the  possibility  of  writing  down  for  the  case  of  interest  the  reactions  that  are  ex¬ 
pressed  by  equations  (3.11).  In  many  cases  it  is  possible  to  write  these  reactions  immediately 
by  intuition.  In  some  cases  it  may  be  necessary  to  formulate  the  reactions  of  equations  (3. 11) by 
applying  to  the  system  under  consideration  the  steps  indicated  by  equations  (3.8)  to  (3.10). 

If  the  number  of  moles  of  the  j-th  component  contained  in  the  hypothetical  system  consisting 
of  components  only  is  denoted  by  q  j,  it  is  evident  that 


ajk  Qk» 
j  =  l 


(3.12) 


k  -  1,  2,  ...  m.  Vie  shall  assume  henceforth  that  the  number  of  components  c  is  equal  to  the  num 
ber  of  different  elements  m.  Introducing  the  ratios  defined  by  equation  (3.1),  we  obtain 


«jk  qj  -  (x^kVx(m))  Qm  =0.  (3.13) 

j=l 
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For  the  systems  to  be  considered  in  this  report,  it  is  convenient  to  assume  that 


qj  -  i» 


(3.14) 


where  the  operator  V  1  denotes  a  summation  over  the  gaseous  components  only.  Equations 


(3.13)  and  (3.14)  constitute  c+1  linear  nonhomogeneous  equations  in  the  quantities  Qm  and  qj, 
j  =  1,  2,  .  .  .  c,  and  they  may  be  employed  to  express  the  qj  in  terms  of  the  ratios  (X^Vx^m^)  and 
thus  in  terms  of  the  parameters  employed  for  the  description  of  the  system.  Equation  (3. 14)  is 
equivalent  to  an  assumption  as  to  the  total  size  of  the  system,  the  weight  MQ  of  the  gaseous 
phase  being  given  by 


qj  mj, 


(3.15) 


where  Mj  is  the  molecular  weight  of  the  j-th  component.  Vt’e  shall  call  the  constants  qj,  deter¬ 
mined  by  means  of  equations  (3.13)  and  (3.14),  the  normalized  stoichiometric  constants. 


Conditions  for  Equilibrium 

The  discussion  of  equilibrium  in  systems  of  many  constituents  assumes  a  particularly  simple 
and  symmetrical  form  in  terms  of  the  definitions  that  have  been  introduced.  U'e  consider  the  gen¬ 
eral  case  where  the  s  substances  are  distributed  among  Tt  coexisting  phases.  The  conservation  of 
mass  in  the  system  as  a  whole  requires  that 


k=l  (  j=l 


c 

X  Hj^Mj 


i-'c+  1 


ni  (k)Mi  [  -  ^  qj  Mj, 

)  r  1 


(3.16) 


where  n^k)  and  nj(k)  are  the  number  of  moles  in  the  k-th  phase  of  the  i-th  and  j-th  substances, 
respectively,  and  Mj  and  Mj  are  the  corresponding  weights  of  1  mole.  As  mass  is  conserved  by 
each  of  reactions  (3.11),  one  may  substitute 


Mj  =  ^  Vij  Mj 

j  =  l 

into  equation  (3. 16).  On  equating  the  coefficients  of  Mj,  there  are  obtained 


(3.17) 


Z.  I  nj (k)  +  z_  vij 

k-1  '  i-c+1 


Vij  nj(k)  =qj, 


(3.18) 
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j  =  1,  2,  .  .  .  c. 

The  derivation  by  means  of  the  method  employed  by  Gibbs, 4/  subject  to  the  constraints  of 
equations  (3.18),  of  the  conditions  for  equilibrium  is  straightforward,  and  it  is  unnecessary  to 
reproduce  it  here.  The  usual  conditions  for  thermal  and  mechanical  equilibrium  are  obtained,  and 
the  conditions  for  chemical  equilibrium  take  the  form, 


^  =  Aj, 


(3. 19) 


M- 


k  =  1,  2,  ...  71; 


-I 

j=l 


Vij  V 


(3.20) 


1,  2,  .  .  .  c;  i  c  +  1,  c  +  2,  .  .  .  s,  where  |i/  ^  and  |i  ^  ^  are  the  chemical 

1  ) 


potentials  in  the  k-th  phase  of  the  i-th  and  j-th  substances,  respectively,  the  A,  j  are  constants 
that  may  be  eliminated  from  equations  (3.19)  and  (3.20),  and  where  it  has  been  assumed  for  sim¬ 
plicity  that  each  substance  of  the  system  is  included  in  every  phase.  If  a  particular  component 
is  excluded  a  priori  from  some  phase,  the  corresponding  member  of  equations  (3.19)  disappears 
as  being  without  meaning,  and  similar  considerations  obtain  if  a  particular  dependent  constituent 
is  excluded  a  priori  from  some  phase. 

Together  with  equation  (3.18),  equations  (3.19)  and  (3.20)  provide,  by  reason  of  their  symme¬ 
try,  the  simplest  basis  for  calculating  the  concentrations  at  equilibrium  of  the  several  constitu¬ 
ents  of  the  system.  In  particular,  if  the  substances  are  limited  to  a  single  phase,  the  equations 
be  come 

s 

vijni=qj,  (3.21) 


nj 


i-c+1 


j=I 


Vij  Uj 


O, 


(3.22) 


j  =  1,  2,  .  .  .  c,  i  =  c  +  l,  ...  s,  for  equilibrium  at  constant  temperature  and  pressure.  Equations 
(3.21)  and  (3.22)  are  also  applicable  to  systems  composed  of  a  gas  phase  in  equilibrium  with  a 
single  pure  solid  phase,  as  in  this  case  it  is  not  necessary  to  retain  explicitly  the  index  label 
of  the  phase. 

The  chemical  potential  of  the  k-th  constituent  of  a  gas  mixture  can  be  expressed  as^ 

o  * 

kk  =  M^k  (T)  RT  log  Pk  •  (3.23) 


4/  Collected  Works  of  J.  Willard  Gibbs:  Longmans,  Green  &  Co.,  New  York  and  London,  1928, 
vol.  1,  pp.  63-82. 

5/  See,  for  example,  Guggenheim,  E.  A.,  Modern  Thermodynamics  by  the  Methods  of  Willard 
Gibbs:  Methuen  &  Co.,  Ltd.,  London,  1933,  P-  72. 


4351 


-  8- 


$  o 

where  p  is  the  fugacity  of  the  k-th  constituent  in  the  mixture  and  |l  (T),  the  chemical  potential 
k  k 


of  the  k-th  constituent  in  the  standard  state  of  unit  fugacity,  is  a  function  of  the  temperature,  T, 
only.  The  chemical  potential  of  a  constituent,  say  the  m-th,  existing  as  a  pure  solid  phase,  can 
be  written  as  a  function  of  temperature  only, 


H  =Ut°  (T),  (3.24) 

m  m 


for  the  range  of  pressures  to  be  considered  in  this  report.  Equations  (3.22)  can  be  written 

*  ^ - -  t  * 

log  Pi  -  /  v  ..  log  p.  =log  Ki(  (3.25) 

where 

c  o  o 

RTlogKi  =  ^  V..  n-  -n-  (3.26) 

j  =1 

Kj,  the  thermodynamic  equilibrium  constant  of  the  reaction  (3- 11)  for  the  formation  of  the  i-th  de¬ 
pendent  constituent  from  the  components,  is  a  function  of  the  temperature  only.  Equations  (3.25) 
are  valid  both  for  homogeneous  gas  mixtures  and  for  heterogeneous  systems  composed  of  a  gas 
mixture  in  equilibrium  with  a  single  pure  solid  phase,  the  summation  being  over  gaseous  com¬ 
ponents  only. 


Degrees  of  Freedom  of  Carbon,  Hydrogen,  Oxygen,  Nitrogen  Systems 

The  phase  rule  of  Gibbs  can  be  obtained  directly  from  Equations  (3*18),  (3*19),  and  (3.20). 
It  can  be  expressed  in  the  form 


f  -  c  +  2  -  n, 

where  f  is  the  number  of  degrees  of  freedom  possessed  by  a  system  of  c  components  and  Tt  phases 
that  depends  on  two  state  variables. 

The  four-component  system  containing  carbon,  hydrogen,  oxygen,  and  nitrogen  consists  of 
two  regions.  There  is  a  region  of  low  oxygen:  carbon  ratio  in  which  the  system  is  heterogene¬ 
ous,  consisting  of  a  gas  phase  in  equilibrium  with  solid  carbon.  If  the  oxygen: carbon  ratio  ex¬ 
ceeds  a  critical  value  (which  depends  on  the  temperature,  pressure,  and  composition),  the  system 
is  homogeneous,  consisting  of  a  single  gas  phase. 


6/  Brinkley,  S.  R.,  Jr.,  Work  cited  in  footnote  3* 
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According  to  the  phase  rule,  the  homogeneous  portion  of  the  four-component  system  pos¬ 
sesses  five  degrees  of  freedom,  and  the  heterogeneous  portion  four  degrees  of  freedom.  The 
state  is  uniquely  determined  by  the  specification  of  two  state  variables,  for  which  purpose  either 
the  temperature  and  pressure  or  the  temperature  and  volume  are  commonly  selected.  The  thermo¬ 
dynamic  requirements  for  the  description  of  the  homogeneous  portion  of  the  system  are  completed 
by  the  specification  of  three  parameters  defining  its  gross  composition.  The  specification  of 
only  two  parameters  defining  gross  composition  is  necessary  in  the  case  of  the  heterogeneous 
portion  of  the  system. 

As  equations  (3.25)  are  independent  of  the  amount  of  solid  phase  actually  contained  by  the 
system,  it  is  convenient  to  define  the  relative  amounts  of  hydrogen,  oxygen,  and  nitrogen  avail¬ 
able  to  the  system.  For  this  purpose,  parameters  D  and  E  of  equations  (3.3)  are  appropriate,  and 
the  specification  of  D  and  E  together  with  two  state  variables  is  sufficient  thermodynamically 
for  the  complete  description  of  the  system.  The  number  of  gram  atoms  of  carbon  present  in  the 
gas  phase  of  the  system  at  equilibrium  (as  carbon  monoxide,  carbon  dioxide,  methane,  etc.)  is 
easily  computed  after  the  equilibrium  composition  is  known.  If  this  quantity  is  denoted  by  Qc, 
we  may  calculate 


Bi=Qh/(4Qc  +  Qh)-  ( 3  -  27) 

As  Qc  is  a  function  of  the  composition  at  equilibrium,  the  ratio  defined  by  equation  (3-27)  is  an 
implicitly  defined  function  of  the  state  (say  the  temperature  and  pressure)  and  of  the  relative 
hydrogen,  oxygen,  and  nitrogen  compositions.  This  functional  relationship  is  expressed  by 


Bf  -3^(0,  E,  p,  T), 


(3.28) 


where  p  is  the  pressure.  Equation  (3.28)  defines  a  surface  in  the  space  with  coordinates  B,  D, 
and  E  that  may  be  termed  the  carbon  boundary  surface.  This  surface  represents  the  limiting  con¬ 
ditions  under  which  a  gas  phase  can  exist  in  equilibrium  with  solid  carbon. 

If  the  relative  gross  composition  of  a  particular  mixture  of  the  four  elements  is  specified  by 
values  of  the  three  quantities  B,  D,  and  E,  as  defined  by  equations  (3.3),  and  if  B  ‘-Bj  fora 
given  temperature  and  pressure,  then  the  mixture  at  equilibrium  consists  of  a  gas  phase  whose 
composition  is  uniquely  determined  by  the  values  of  D  and  E  in  equilibrium  with  Qc  -  Qc  gram 
atoms  of  carbon,  where 


(3.29) 


As  the  ratios  D  and  E  can  be  expressed  as  functions  of  A,  B,  and  C  by  means  of  equations 
(3-4) 

D  =  D  (A,  B,  C), 

E  =  E  (B,  C), 
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it  follows  that  the  functional  relation,  equation  (3.28),  defining  the  carbon  boundary  surface,  can 
be  converted  to  the  form, 

Bj  =Bi(A,  C,p,T),  (3.30) 

and  inverted  to  the  form 

Cl  =CX  (A,  B,  p,  T)  (3.31) 

by  suitable  numerical  procedures.  Here  Cj  is  the  critical  value  of  C  for  the  carbon  boundary  sur¬ 
face,  the  heterogeneous  portion  of  the  system  corresponding  to  the  region  C  <(21.  The  represen¬ 
tation,  equation  (3.28),  is  the  natural  one  to  employ  for  investigations  within  the  heterogeneous 
portion  of  the  system,  and  the  representation,  equation  (3.31),  is  the  more  convenient  one  to  em¬ 
ploy  for  the  description  of  the  extent  of  the  homogeneous  region.  We  shall  show  in  a  later  sec¬ 
tion  that  it  is  possible  to  determine  the  relation  (3-31)  by  direct  calculation  for  given  values  of 

A,  B,  p,  and  T,  as  an  alternative  procedure  to  the  series  of  interpolations  required  to  obtain  the 
representation  (3*31)  from  the  representation  (3.28). 

If  C  >Cj  (A,  B,  p,  T),  the  system  consists  of  an  homogeneous  gas  phase  with  five  degrees 
of  freedom.  The  complete  description  of  the  system  therefore  requires  the  specification  of  two 
state  variables  and  three  parameters  describing  the  gross  composition  of  the  system,  the  set  A, 

B,  C  of  equations  (3-2)  being  a  particularly  convenient  set  for  the  consideration  of  combustion 
gases. 

IV.  CALCULATION  OF  THE  EQUILIBRIUM  COMPOSITION 


The  composition  of  a  system  at  chemical  equilibrium  is  easily  calculated  where  there  is  only 
a  single  reaction  (one  dependent  constituent)  to  be  considered.  In  this  case,  the  concentration 
of  each  constituent  can  be  related  to  a  single  variable,  "the  degree  of  reaction,”  and  the  solu¬ 
tion  of  the  mass  action  equation  is  straightforward.  Difficulties  are  encountered  if  this  method 
is  extended  to  a  consideration  of  two  simultaneous  equilibria;  and  when  the  number  of  such  si¬ 
multaneous  equilibria  becomes  large,  the  ordinary  methods  become  very  laborious.  In  a  recent 
publication,?/  a  systematic  procedure  for  calculating  the  equilibrium  composition  of  a  system  of 
many  constituents  was  presented.  This  method  presents  a  simple  rule  for  formulating  the  work 
program  of  such  calculations,  with  the  result  that  very  little  time  is  required  for  setting  up  a  par¬ 
ticular  problem.  The  systematic  nature  of  the  procedure  makes  the  method  well  adapted  to  com¬ 
putation  by  modern  automatic  computational  equipment.  In  the  publication  cited,  the  method  was 


l]  Brinkley,  S.  R.,  Jr.,  Calculation  of  the  Equilibrium  Composition  of  Systems  of  Many  Constit¬ 
uents:  Jour,  Chem.  Phys.,  vol.  15,  1947,  pp.  107-110. 
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developed  for  systems  of  a  very  general  nature.  In  subsequent  communications, — /  applica¬ 
tion  of  the  method  was  restricted  to  the  calculation  of  the  equilibrium  composition  of  mixtures 
consisting  of  a  single  homogeneous  gas  phase  or  of  a  gas  phase  in  equilibrium  with  a  single  pure 
solid  phase,  and  it  was  assumed  that  the  gas  phase  is  adequately  described  by  the  ideal  gas 
equation  of  state.  By  taking  advantage  of  these  restrictions,  it  was  possible  to  formulate  a  com¬ 
putational  method  applicable  to  these  particular  cases,  which  is  substantially  simpler  and  more 
systematic  than  the  more  general  method.  In  addition,  it  was  possible  to  formulate  a  special 
method,  applicable  in.  restricted  cases,  that  has  the  advantage  of  converging  considerably  more 
rapidly  than  the  more  systematic  general  method. 

In  order  that  the  present  discussion  may  be  self  contained,  we  repeat  portions  of  the  develop¬ 
ment  of  the  basic  equations  that  have  been  published  elsewhere. 


The  General  Method 


The  formal  methods  of  the  preceding  section  have  resulted  in  a  statement  of  the  conditions 
for  equilibrium  in  a  form  that  is  particularly  well  suited  to  serve  as  a  basis  for  the  computation 
of  the  equilibrium  composition.  We  shall  henceforth  limit  consideration  to  the  cases  of  a  homo¬ 
geneous  gaseous  system  or  of  a  heterogeneous  system  composed  of  a  gas  mixture  in  equilibrium 
with  a  single,  pure,  solid  phase. 

The  conservation  of  each  element  was  shown  to  require  that 

s 

nj  +  /*  v  ij  ni  =  qj>  (4- *) 

i=c+  1 

j  =  1,  2,  .  .  .  c,  where  nj  and  rij  are  the  number  of  moles  of  the  i-th  dependent  constituent  and 
j-th  component,  respectively,  in  the  equilibrium  mixture  containing  Mc  grams  of  gas,  qj  is  the 
number  of  moles  of  the  j-th  component  in  the  hypothetical  mixture  containing  components  only 
and  where  V  -  is  the  coefficient  to  the  symbol  of  the  j-th  component  in  the  equation  for  the  for¬ 
mation  of  the  i-th  dependent  constituent  (involving  only  components  as  reactants).  The  con¬ 
stants  qj  are  normalized  in  accord  with  the  relation, 


8/  Kandiner,  H.  J.,  and  Brinkley,  S.  R.,  Jr.,  Calculation  of  Complex  Equilibrium  Relations:  Ind. 
Eng.  Chem.,  vol.  42,  May  1950,  pp.  850-855. 


9/  Brinkley,  S.  R.,  Jr.,  and  Smith,  R.  W.,  Jr.,  Calculation  of  the  Equilibrium  Composition  of 
Systems  of  Many  Constituents:  Proc.  Scientific  Computation  Forum  1948,  International 
Business  Machines  Corp.,  New  York,  1950,  pp.  77-82. 

10/  Brinkley,  S.  R.,  Jr.,  and  Smith,  R.  W.,  Jr.,  Calculation  of  the  Equilibrium  Composition  of 
Homogeneous  Multicomponent  Systems:  Proc.  Seminar  on  Scientific  Computation  1949, 
International  Business  Machines  Corp.,  New  Tork,  1950,  pp.  58-63- 
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where  the  operator  1  denotes  a  summation  over  gaseous  components  only,  and  this  normal- 

-i 


ization  implies  that  the  weight  of  the  gas  phase  M0  is  assigned  the  arbitrarily  assumed  value 
given  by  equation  (3.15). 


The  conditions  for  chemical  equilibrium,  given  in  equations  (3.25),  can  be  put  in  the  form, 


(4.3) 


i  =c+l,  .  .  .  s,  where  pj  and  pj  are  the  fugacities  of  the  i-th  dependent  constituent  and  j-th  com¬ 
ponent,  respectively,  Kj  is  the  thermodynamic  equilibrium  constant,  a  function  of  the  temperature 
only  that  is  defined  by  equation  (3.26),  of  the  reaction  for  the  formation  of  the  i-th  dependent  con¬ 
stituent,  and  where  the  operator  7t  j  denotes  that  the  repeated  product  is  to  be  taken  over  gaseous 
components  only.  The  equilibrium  constant  Kj  is  most  easily  obtained  by  the  relation, 


(4.4) 


where  Kij  and  Kfj  are  the  equilibrium  constants  for  the  formation  of  the  i-th  dependent  constituent 
and  j-th  component,  respectively,  from  the  elements. 

For  each  gaseous  component  we  employ  equations  (4.1)  in  the  form 


X:  + 


XT 

i=C+  l 


Vij  xi 


qj/n, 


(4.5) 


where  Xj  and  Xj  are  the  mole  fractions  in  the  equilibrium  mixture  of  the  i-th  dependent  constituent 
and  the  j-th  gaseous  component,  respectively,  and  n  is  the  total  number  of  moles  of  gas  in  the 
equilibrium  mixture  corresponding  to  the  normalized  constants  qj.  In  view  of  equation  (4.2) 


where 


V  j  xj  -  1/n, 


(4.6) 


The  mole  fractions  of  the  gaseous  constituents  are  subject  to  the  identity  relation 
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Therefore,  equation  (4.6)  becomes 


s 

l/n  =  1  +  T>~  (  V  i-1)  xi, 

i=t+  1 


and  equations  (4.5)  may  be  written 


Xi  =C>j 


i~c+  1 


[ vij  -  qj  (vi  -  i)  1  xi. 


(4.7) 


there  being  a  member  of  the  set  (4.7)  for  each  gaseous  component. 

If  an  activity  coefficient  fk  of  the  k-th  gaseous  constituent  is  defined  by  the  relations 

Pk  =xk  fk  P> 


Lim  (p->0)  fk  =  1, 

equations  (4.3)  can  be  written  in  the  form 

Xi  =ki  nj  Xj  V‘j 

i  =  c  +  1,  .  .  .  s,  where  the  ki  are  functions  of  temperature  and  pressure  given  by 

Vj-1 

ki  “Ki  P  ri. 


(4.8) 


(4.9) 


(4.10) 


in  which  the  Tj  are  in  general  functions  of  the  temperature  and  pressure.  In  this  section  we  shall 
assume  that  the  gas  mixture  is  ideal  and  therefore  that  Pj  =  1  for  every  i.  In  a  later  section  we 
shall  consider  the  evaluation  of  the  Pj  for  imperfect  gas  mixtures. 

The  computation  of  the  equilibrium  composition  requires  the  simultaneous  solution  of  equa¬ 
tions  (4.7)  and  (4.9.  If  xj«Xj  for  all  i  and  j,  the  solution  may  be  carried  out  by  a  simple  itera¬ 
tion  method.  U/  An  approximate  set  of  values  is  chosen  for  the  Xj.  (In  the  absence  of  any  cri¬ 
teria  for  the  choice  of  the  initial  set,  one  may  take  Xj  =q j.)  Equations  (4.9)  are  employed  in 
the  computation  of  corresponding  values  of  the  Xj.  These,  in  turn,  are  employed  with  equations 
(4.7)  for  the  determination  of  an  improved  set  of  values  for  the  Xj.  This  iterative  process  is  con¬ 
tinued  until  the  difference  between  successive  approximations  to  the  Xj  is  less  than  the  desired 
precision  of  the  computation. 

1 1/  Scarborough,  J.  B.,  Numerical  Mathematical  Analysis:  Johns  Hopkins  Press,  Baltimore,  Md., 
1930,  pp.  191-195. 
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The  convergence  of  this  simple  iteration  method  is  very  slow  for  larger  relative  values  of  the 
xj,  and  when  the  xj  and  Xj  are  of  the  same  order  of  magnitude,  this  method  may  not  converge  at  all 
for  any  choice  of  components.  A  more  powerful  computational  procedure  is  provided  by  the 
Newton-Raphson  method.^/  Equations  (4.7)  may  be  written  in  the  form 

Fj  =0  (4-  n) 

for  each  gaseous  component,  where 

s 

Fj  =  qj  G  '  *j  ‘ 

i=c+  1 


s 


i=c+ 1 


1)  Xj. 


We  seek  the  solution  of  equations  (4.11)  subject  to  equations  (4.9).  If  the  functions  E:  are  ex¬ 
panded  in  Taylor  series  about  an  approximate  set  of  values  oi  the  variables  Xj  with  neglect  of 
terms  involving  derivatives  of  second  and  higher  orders,  there  results  a  set  of  linear  equations 
that  can  be  compactly  represented  in  the  notation  of  matrices  by 

Pn’] 

where  the  r-th  and  (r  +  I)-th  approximations  to  the  composition  are  related  by 


(r) 


(4.12) 


and  where  the  elements  of  the  matrix  are  given  by 

V  V 

where 


U)i'  =Xj  6)l 


s 

>'  +  21  vij  Vi)  xi’ 


i-c+  1 


(4. 13) 


(4. 14) 


s 


i  -c+  1 


and  where  Ujj  >  —  Uj  ' j .  The  superscript  r  indicates  that  the  designated  quantity  is  to  be  evaluated 
with  the  r-th  approximation  to  the  composition  of  the  system,  and  6jj  t  is  the  Kronecker  delta. 


12/  Scarborough,  J.  B.,  Work  cited  in  footnote  11,  pp.  178  and  187. 
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Criteria  for  the  choice  of  components  that  results  in  the  most  rapid  convergence  of  the  itera¬ 
tion  process  can  be  developed  from  the  remainders  to  the  two-term  Taylor  series  expansions  of 
functions  Fj.  However,  the  resulting  expressions  are  too  cumbersome  for  practical  utility,  and 
in  practice  the  convergence  will  be  found  to  be  satisfactory  if  the  components  are  selected  so  as 
to  minimize  the  quantities  kj,  i  =  c  +  1,  .  .  .  ,  s. 

A  Special  Method 


An  alternative  computational  procedure  is  very  useful  in  certain  special  cases.  We  employ 
equations  (4.  1), 


nj  +  V  ij  n;  -  qj, 

i=c+  1 


(4. 1) 


and  write  equations  (4.9)  in  the  form 


T  V: 


nj  -kj  (n)  rt:  n:  , 


(4.15) 


ki  (n)  -kj  n 


1-  vi 


i  =c  +  1,  ...  s,  by  means  of  the  substitution  n^  -xj^n.  In  view  of  equations  (4.  1)  and  (4.2),  the 
total  number  of  moles  of  gas  n  is  given  by 


n  =  1  +  (  1  -  Vj)  ni. 

i=t+  1 


(4. 16) 


We  define  a  variable  |i  by 


'  r'j 

li  -  n  j  nj  , 


(4. 17) 


and  rewrite  equations  (4. 15)  in  the  form 

nj  =kj  (n)M.  "  TZ-  Rj 


.  £j  -  Gjj 


(4.  18) 


°ij  =  vij  -  nj. 

i  =c+  1,  .  .  .  s.  We  endeavor  to  choose  the  T\j,  so  that  tor  each  i  there  exists  only  a  single  value 
of  j,  say  j  ',  for  which  G  jj  does  not  vanish,  and  so  that  for  some  value  of  j,  say  j*,  Gjj*  vanishes 
for  every  value  of  i.  If  this  can  be  done,  equations  (4.18)  reduce  to  expressions  of  the  form 


i  <$j  ®ij  ' 

nj  =kj  (n)|i  nj  i  , 


(4. 19) 


i  -  c  +  1,  .  .  .  s,  j  -j*. 
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Introduction  of  the  variable  d  leads  immediately  to  a  separation  of  variables  in  the  set  of 
equations  composed  of  equations  (4.19)  and  those  of  equations  (4.1),  for  which  / j*  when  mass 
action  expressions  of  the  form  (4. 19)  can  be  obtained  with 

T) :  * 

U=n-*  '  (4.20) 


Then  T^*  is  the  lowest  common  divisor  of  those  coefficients  Vjj*  that  are  nonzero  and  Ojj  =  Vjj 
for  j  7^  j*.  By  substitution  of  equations  (4. 19)  into  equations  (4. 1),  we  obtain 


'  k  vii 

vij  kj  (n)|i  nj  -  9j  =  O 


(4.21) 


for  j  f  j*,  which  can,  at  least  in  principle,  be  solved  for  nj  by  ordinary  algebraic  methods  for 
given  values  of  n  and  |i.  If  equations  (4.19)  and  (4.20)  are  then  employed,  a  solution  can  be  found 
in  the  form 


nj  -  nj  (n,  |l), 
ni  =  nj  (n,  d), 


j  -  1,  2,  .  .  .  c;  i  =  c  +  1,  .  .  .  s. 

The  values  of  n  and  |i  leading  to  an  exact  solution  of  the  equilibrium  conditions  must  satisfy 
equation  (4.16)  and  the  member  of  equations  (4.1)  for  which  j  =  j*.  These  can  be  written  in  the 
form 


F(U)=0,  F(n)=G, 


(4.22) 


F  (d)  =qj* 


s 

'  nj*  ■  ^  vij  *  ni’ 


i  -c+ 1 


F  (n)  -n  -  1 


2 

i-c+  1 


(Vj  -  1)  nj 


If  the  functions  F  (d)  and  F  (n)  are  expanded  in  Taylor  series  about  approximate  values  of  n  and 
0  with  neglect  of  terms  involving  second  and  higher  derivatives,  there  result  a  pair  of  linear 
equations, 


S’  N. 

(r) 

A11 

(r) 

A12 

(D 

h(d)- 

= 

(r) 

F  (d) 

(r) 

(r) 

(f) 

(r) 

^21 

a22^ 

h  (n) 

b«.  / 

= 

1  F  (n)  , 

2  2 

where  the  r-th  and  (r+  l)-th  approximation  to  d  and  n  are  related  by 
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(4.24) 


(r+1) 

(r) 

[i  + 

(0, 

=  M. 

H  (M.)  ], 

(r+1) 

(r) 

(r) 

-  n 

[i  + 

h (n)  ]  , 

and  where  the  coefficients  are  given  by 


A 


11 


L  + 


i-c+  1 


vij*  ^i  ni 


vij*  vij(i)  Uj(i)  ni> 


l12 


i-c+  1 


vij*  "  vi)  ni  +  /  vij*  vij(i)  ^j(i)  ni’ 


i=c+  l 


(4.25) 


A21  ~nj  M  +  ^  (1  -  Vj)  fii  ^  (1  -  Vj)  V ij(£)  Uj(i)  q, 

i  =c  + 1  i  =c  + 1 


s 


i=c+  1 


(1  -  Vi)2 


S 

+  (1-  VA)  Vij(i)  Vj(i)  ni,. 

i=c  +  l 


where  j(i)  is  the  value  of  j,  for  which  O  ij  does  not  vanish  for  a  given  i,  and  where 

L  =  1/n-  M  =0, 


1 

II 

1  Vij  £i  "i 

nj  (Vij)2ni 

_i(j) 

L  Kj)  _ 

vi =- 

2  vij(i*vi)  ni 

nj  +  (vij)2ni 

K]) 

«j) 

(4.26) 


The  superscript  r  indicates  that  the  designated  quantity  is  to  be  evaluated  with  the  r-th  approxi¬ 
mation  to  the  composition  of  the  system. 

Introduction  of  the  variable  (_l  also  leads  to  a  separation  of  variables  in  the  set  of  equations 
composed  of  equations  (4.  19)  and  those  of  equations  (4.  1),  for  which  j  f  j*  when  mass  action  ex¬ 
pressions  of  the  form  (4.19)  can  be  obtained  with 


M-  =nj  t  71  i'  nj*,  (4.27) 

j'  fj*.  Then  £j=Vjj*,  Ojj  =  V  jj  for  j  ^  j j*,  and  C  jj  t  =  Vjj »  -  V^j*  T\j  i.  The  stoichiometric 
conditions  for  j  f)*  can  then  be  put  into  a  form  analogous  to  equations  (4.21) 
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^  '  £.  v 

nj  +  Z.  Vjjkj^p1^  li  -  qj  =  0,  j  ii\ 


1  vl) 

i(j) 


-  nj'  V" 

nj  i-  rij  f  M-  nj  i  + 


i(j  ') 


t  °ii ! 

0  jj  tki  (n)M.  nj  -qji  +  Tlji  qj*  =0, 


(4.28) 


which  can,  at  least  in  principle,  be  solved  for  nj,  j  f1  j* ,  by  ordinary  algebraic  methods  for  given 
values  of  n  and  pi.  If  equations  (4.19)  and  (4.27)  are  then  employed,  a  solution  can  be  found  in 
the  form 

nj  =  nj  (n,  p), 
nj  =  nj  (n,  |i), 

j  =  1,  2,  .  .  .  c;  i  =  c  +  1,  .  .  .  s. 

The  values  of  n  and  p  leading  to  an  exact  solution  of  the  equilibrium  conditions  must  satisfy 
equations  (4.22),  The  application  of  the  Newton-Raphson  method  to  obtain  an  improved  pair  of 
values  of  n  and  p  requires  the  solution  of  equations  (4.23)  and  the  application  of  equations  (4.24). 
The  coefficients  of  equations  (4.23)  are  given  by  equations  (4.25)  but  with 


L  =  1  -  r\j  i  Uj  t  ,  M  =  -  i3j  i  Vj  i 


1 

II 

pr 

2  vij  h  ni 

nj  +  2  <vij>2  ni 

Li(j)  J 

i(j)  J 

Uj’  = 

—  _  T 

r'j  1  nj*  -  a  jj  T  <5j  nj 

nj  i  +  Tlj'nj*  +  (Ojji)2nj 

L  i(j')  J 

i(j  ') 

v- 

X  v  jj  (1  -  V  j)  nj 

nj  +  X  ( v ij) 2  ni 

_J(j) 

Ki) 

Ojji  (1  -  Vj)  nj 


2 

nj  i  +  r\ j  nj  * 


+ 


(°ij)2 


replacing  equations  (4.26). 


-1 


(4.29) 


The  special  method  converges  considerably  more  rapidly  than  the  general  method  and  can  be 
advantageously  employed  for  computations  where  it  is  applicable.  It  has  the  disadvantage,  how¬ 
ever,  that  the  computation  program  for  a  single  iteration  is  longer  and  less  systematic  than  for 
the  general  method.  The  final  results  are  easily  expressed  in  terms  of  mole  fractions  by  the 
application  of  the  relation 

xk  =nk/n.  (4.30) 
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Modifications  for  Computation  at  a  Specified  Density 


Many  of  the  applications  of  thermodynamics  to  combustion  processes  are  such  that  it  is  nat¬ 
ural  to  specify  the  state  by  means  of  the  variables  temperature  and  pressure.  The  methods  for 
the  calculation  of  the  equilibrium  composition  that  have  been  described  assume  that  the  temper¬ 
ature  and  pressure  have  been  explicitly  specified,  i.e.,  that  the  equation  of  state  i s  of  the  form 

p/  p  =nRT  cp  (p,T)/M0, 

where  n  is  the  number  of  moles  of  gas  contained  in  MQ  grams  of  gas,  P  is  the  density,  and  *P  (p,T) 
is  explicit  in  the  pressure  and  temperature  and  unity  for  the  ideal  gas. 

In  certain  applications,  as  in  those  of  internal  ballistics,  it  is  more  natural  to  specify  the 
state  by  means  of  the  variables  temperature  and  density  (or  volume).  For  these  cases  the  equa¬ 
tion  of  state  may  be  assumed  to  be  of  the  form 

p/p  =nRT  \|f  (p,T)/M0, 

where  (  P,T)  is  explicit  in  the  density  and  temperature  and  unity  for  the  ideal  gas.  The  meth¬ 
ods  for  the  calculation  of  the  equilibrium  composition  at  a  specified  density  are  analogous  for 
those  at  a  specified  pressure,  but  they  differ  in  certain  details. 

The  conditions  for  equilibrium,  equations  (4.3),  can  be  advantageously  written  in  the  form, 


1  !  V:: 

^  =  kj  Tlj  nj  ), 


i  =c+l,  .  .  .  s,  where  the  k  j  are  functions  of  temperature  and  density  given  by 


(4.31) 


k|  =  Kj  (rtp/Mq)  Vi'1  rj, 

1  (4.32) 

rj  =  [^(P,T)]  Vi'lfi’1rtj  fj  Vij, 

in  which  the  F  j  are  in  general  functions  of  the  temperature  and  density  and  are  unity  for  every  i 

if  the  gas  mixture  is  ideal.  In  a  formal  way,  equations  (4.31)  and  (4.32)  can  be  obtained  from  the 

analogous  expressions,  equations  (4.9)  and  (4.10),  employed  in  the  general  method  above  by  re- 

1  ! 

placing  Xj  and  xj  by  nj  and  nj,  respectively,  and  by  replacing  kj  by  kj.  We  note  that  the  kj  of 
equations  (4.15)  and  (4.  31)  are  identical,  although  they  are  evaluated  in  different  ways. 

Equations  (4.31)  determine  the  concentrations  of  the  dependent  constituents  for  assumed 
values  of  the  concentrations  of  the  components.  In  order  that  these  calculations  lead  to  an  exact 
calculation  of  the  equilibrium  composition,  it  is  necessary  that  the  results  satisfy  the  stoichio¬ 
metric  conditions,  equations  (4.1),  which  we  write  in  the  form 

F  j  =0,  (4.33) 
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where 


s 


i-c+  1 


In  a  formal  way,  equations  (4.33)  can  be  obtained  from  the  analogous  expressions,  equations 
(4. 11)  by  replacing  Xj  and  xj  by  nj  and  nj,  respectively,  and  by  setting  G  =  1. 


The  application  of  the  Newton-Raphson  method  to  the  improvement  of  an  approximate  set  of 
values  of  the  concentrations  of  the  components  requires  the  solution  of  equations  analogous  to 
equations  (4. 12), 


s'  \ 

s'  \ 

f'  N 

,(r) 

,(r) 

l  (f) 

A  jj' 

h  j- 

F  j 

\  > 

\  X1 

where  the  r-th  and  the  (r+ l)-th  approximations  to  the  composition  are  related  by 


nj(r+1>  =nj(r)(l  +  h'(t)); 


(4.35) 


and  where  the  elements  of  the  matrix  are  given  by 


A'jj'  =nj  v  +2l  vij  v  ni- 


i-c+  1 


(4.36) 


In  a  formal  way,  equations  (4.34),  (4.35),  and  (4.36)  can  be  obtained  from  the  analogous  expres¬ 
sions,  equations  (4.12),  (4.13),  and  (4.14),  respectively,  by  replacing  Xj  and  Xj  by  nj  and  nj,  re¬ 
spectively,  and  by  setting  Vj  =0  for  all  j. 

It  is  evident  that  the  general  method  for  computations  at  specified  density  is  similar  in  out¬ 
line  with  the  previously  described  general  method  for  computations  at  specified  pressure,  and 
that  a  particular  computational  program  for  the  former  case  can  easily  be  derived  in  a  formal  way 
from  the  equivalent  program  for  the  latter  case.  The  total  number  of  moles  n  contained  in \10 
grams  of  gas  mixture  can  be  calculated  with  equation.  (4.  16). 

The  alternative  computational  method,  which  we  have  called  the  special  method,  is  easily 
applied  without  substantial  modification  to  calculations  at  specified  density.  Vie  again  employ 
the  stoichiometric  conditions  as  expressed  by  equations  (4.1)  and  the  equilibrium  conditions  as 
expressed  by  equations  (4.15).  The  constants  k 'j  appearing  in  the  equilibrium  conditions  are  to 
be  calculated  by  equations  (4.32),  and  tor  specified  temperature  and  density  they  are  independent 
of  the  total  number  of  moles  n.  Consequently,  it  is  possible  to  seek  a  solution  of  the  form 

nj  =nj(u) 
nj  =nj(|i) 

j  —  1,  2,  .  .  .  c;  i  =  c+  1,  .  .  .  s,  for  those  cases  where  the  introduction  of  a  variable  u  results  in  a 
separation  of  variables. 
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The  value  of  |i  leading  to  an  exact  solution  of  the  equilibrium  conditions  must  satisfy  the 
condition, 


t(U)  =o 


(4.37) 


where  F(d)  is  defined  by  equations  (4.22).  The  application  of  the  Newton-Raphson  method  to 
obtain  an  improved  value  of  [!  requires  the  application  of  the  relation 

|i(r+  1)  =  [j_( r)  [  1  +  h  (M.)  <r>  (4.38) 

where 

(r)  (r)  (r) 

h(u)=F(u)  /An  (4.39) 

and  where.  Ajj  has  been  defined  as  one  of  equations  (4.25).  Equations  (4.37),  (4.38),  and  (4.39) 
are  obtained  from  the  analogous  expressions,  equations  (4.23),  (4.24),  and  (4.25),  respectively, 
by  retaining  only  those  terms  of  significance  in  the  latter  equations  after  setting  F(n)  =0,  h(n)  =0. 


Calculation  of  the  Equilibrium 


Composition  of  Imperfect  Gas  Mixtures 


In  the  previous  discussion  it  has  been  assumed  that  the  gas  mixture  can  be  described  ade¬ 
quately  by  the  ideal  gas  equation  of  state.  At  the  temperatures  of  interest  in  the  consideration 
of  combustion  gases  this  assumption  is  adequate  for  most  purposes.  However,  for  combustion 
processes  occurring  at  high  pressures,  as  in  closed  bombs  or  in  gun  barrels,  the  ideal  gas  equa¬ 
tion  of  state  is  no  longer  adequate,  and  it  is  necessary  to  consider  the  effects  of  gas  imperfec¬ 
tion  on  the  calculation  of  the  equilibrium  composition.  In  the  present  section,  we  indicate  the 
modifications  to  the  previously  described  computation  methods  that  are  necessary  to  take  account 
of  gas  imperfections,  and  we  assemble  the  relations  required  for  this  purpose.  A  comprehensive 
discussion  of  the  evaluation  of  the  thermodynamic  properties  of  mixtures  of  real  gases,  with 
derivations  of  the  necessary  relations,  has  been  given  by  Beattie.  15/ 

For  calculations  that  are  carried  out  at  a  specified  pressure,  the  mass  action  constant  kj  is 
related  to  the  thermodynamic  equilibrium  constant  Kj  by  equation  (4.10), 

ki  =Ki  P  Jl  1  ^i’ 

where 

r=  =f:*1  71  '  f  V  ij,  (-4-40) 

j  ) 

and  where  the  activity  coefficient  of  the  k-th  constituent  is  the  ratio  of  the  fugacity  to  the  partial 
pressure  of  the  k-th  constituent  in  the  gas  mixture 

fk  =Pk/xkP, 


13/  Beattie,  J.  A.,  The  Computation  of  the  Thermodynamic  Properties  of  Real  Gases  and  Mix¬ 
tures  of  Real  Gases:  Chem.  Rev.,  vol.  44,  February  1949,  pp.  141-192. 


tending  to  the  unity  as  the  pressure  tends  to  zero.  In  the  general  case  of  the  non-ideal  gas,  the 
Tj  are  functions  of  temperature,  pressure,  and  composition. 


If  the  equation  of  state  is  explicit  in  the  pressure  and  of  the  form 

P  _  nRT 


<P(p,T)> 


(4.41) 


P  »'o 

where  lP.(p,T)  is  generally  a  function  of  the  composition  as  well  as  of  the  temperature  and  pres¬ 
sure,  the  activity  coefficient  of  the  k-th  constituent  of  the  gas  mixture  is  given  by 


log  f.  ~J  [r-P-  1  +  n 


d  qp 


on, 


1  dp 
J  P  ' 


(4.42) 


This  expression  can  be  evaluated  explicitly  for  a  particular  equation  of  state  and  employed  to 
determine  the  activity  coefficients  of  the  constituents  of  the  mixture  as  functions  of  temperature, 
pressure,  and,  in  general,  composition.  Equation  (4.40)  can  then  be  employed  to  evaluate  the 
functions  Pj  that  determine  the  effect  of  gas  imperfection  on  the  mass  action  constants  kj. 

If  the  activity  coefficients  are  functions  of  the  composition  of  the  system,  the  Pj  must  be 
evaluated  with  each  successive  approximation  to  the  composition.  The  Pj  are  relatively  insensi¬ 
tive  to  the  composition,  and,  consequently,  it  is  unnecessary  to  modify  the  computation  proce¬ 
dure  by  means  of  which  a  particular  approximation  to  the  composition  is  determined  by  the  one 
preceding  it. 

Lewis  and  Randall!^/  have  proposed  as  a  basis  for  an  approximate  thermodynamic  treatment 
of  imperfect  gas  mixtures  the  rule 

Pk=xkpk(p)>  (4.43) 


where  Pk  (p)  is  the  fugacity  of  k-th  constituent  as  a  pure  gas  at  the  temperature  and  total  pres¬ 
sure  of  the  mixture.  If  thre  k-th  constituent  obeys  the  gas  law 

pv  =nkRT  cpk(p,T),  (444) 


where  9  k  is  a  function  of  pressure  and  temperature  only  and  V  is  the  volume  occupied  by  nk 
moles  at  temperature  T  and  pressure  p,  the  Lewis  and  Randall  rule  results  in  the  expression 

p  dp 

lo8  *k  =  E<Pk  '  "p~>  (4.45) 

for  the  activity  coefficient  of  the  k-th  constituent.  Where  applicable,  equation  (4.45)  is  more  con¬ 
venient  than  the  exact  form,  equation  (4.42),  as  it  is  independent  of  the  composition  of  the  mix¬ 
ture.  The  range  of  validity  of  the  approximate  expression  for  a  particular  application  may  be 
determined  by  comparison  between  results  obtained  from  equations  (4.42)  and  (4.45)  with  an  ex¬ 
plicit  form  of  the  equation  of  state. 


14/  Lewis,  G.  N.,  and  Randall,  M.,  Thermodynamics:  McGraw-Hill  Book  Co.,  Inc.,  New  York, 
1923,  pp.  225-227. 


For  calculations  at  a  specified  density,  the  equation  of  state  may  be  assumed  to  be  of  the 
form 


£ 

P 


nKT 


V  (P,T), 


(4.46) 


where  W  (  P,T)  is  generally  a  function  of  composition  as  well  as  of  temperature  and  density.  The 
mass  action  constants  kq  for  this  case  have  been  related  to  the  corresponding  thermodynamic 
equilibrium  constants  by  the  expressions 


wi  th 


ki  =Ki 


RTP 

M, 


V  :- 1 


r 


T  V :  - 1  .  V::  ,  V:: 

T  i  =  [ty(P,T)]  fi'lrCjfj  =8i'lrtjgj 
where  we  have  introduced  an  alternative  activity  coefficient, 


gk  =  fk  T  (P,T). 


(4.47) 


(4.48) 


The  activity  coefficient  gk  of  the  k-th  constituent  of  the  mixture  is  given  by  the  general  thermo¬ 
dynamic  relation 


,p  ot  d  o 

log  gk  =  i  [v-l  +  ngr]  (4.49) 

O  K  p 

This  expression  can  be  evaluated  explicitly  for  a  particular  equation  of  state  and  employed  to 
determine  the  activity  coefficients  of  the  constituents  of  the  mixture  as  functions  of  temperature, 
density,  and,  in  general,  composition.  Equations  (4.47)  can  then  be  employed  to  evaluate  the 
functions  Pj  that  determine  the  effect  of  gas  imperfection  on  the  mass  action  constants  kj  . 

n  t 

If  the  activity  coefficients  are  functions  of  the  composition  of  the  system,  the  i  j  must  be 
evaluated  with  each  successive  approximation  to  the  composition.  The  1/  are  relatively  insen¬ 
sitive  to  the  composition,  and  it  is  unnecessary  to  modify  the  computation  procedure  by  means  of 
which  a  particular  approximation  to  the  composition  is  determined  by  the  one  preceding  it. 

As  a  basis  for  the  approximate  thermodynamic  treatment  of  gas  mixtures,  GibbslS/  has  pro¬ 
posed  that  the  pressure  of  a  mixture  of  gases  be  assumed  to  be  equal  to  the  sum  of  the  pressures 
of  the  different  gases  as  existing  each  by  itself  at  the  same  temperature  and  with  the  same  value 
of  its  chemical  potential.  If  the  k-th  constituent  obeys  the  gas  law 

pvk  =RT  ¥k(T,vk),  (4.50) 

where  vk  is  the  molar  volume  of  the  k-th  constituent,  then  the  assumption  of  Gibbs  has  been 
shown  by  Beattie!^/  to  lead  to  the  approximate  relation 


15/  Gibbs,  J.  W.,  Collected  Works:  Longmans,  Green  &  Co.,  New  York,  1906,  vol.  I,  pp.  155-158. 
16/  Beattie,  J.  A.,  Work  cited  in  footnote  13. 
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(4.51) 


r  P  r  n  dp 

log  gk  "  /  t¥k  ‘  1  +  «k  'dnk  ]  “ 

for  the  activity  coefficient  of  the  k-th  constituent.  This  expression  is  not  in  general  independent 
of  the  composition  of  the  mixture.  It  may,  however,  lead  to  simpler  expressions  than  those  de¬ 
rived  from  the  general  relation,  equation  (4.49).  The  range  of  validity  of  the  approximate  expres¬ 
sion  for  a  particular  application  may  be  determined  by  comparison  between  results  obtained  from 
equations  (4.49)  and  (4.5  1)  with  an  explicit  form  of  the  equation  of  state. 

V.  COMPUTATION  PROGRAMS  FOR  CARBON-HYDROGEN-OXYGEN-NITROGEN  SYSTEMS 


A  large  number  of  the  fuel-oxidant  systems  of  scientific  and  technical  interest  yield  on  com¬ 
bustion  a  mixture  of  products  composed  of  compounds  of  carbon,  hydrogen,  oxygen,  and  nitrogen. 
Because  of  the  importance  of  such  systems,  we  give  explicitly  in  this  section  the  relations  ob¬ 
tained  by  the  application  of  the  general  methods  of  the  preceding  section.  These  computational 
programs  will  also  serve  to  illustrate  the  application  of  the  general  methods  to  the  calculation 
of  the  equilibrium  composition  of  particular  systems. 

The  computation  programs  will  be  written  for  calculations  at  specified  pressure.  We  do  not 
present  explicidy  the  relations  for  calculations  at  specified  density.  We  have  observed  in  a  pre¬ 
ceding  section  that  the  relations  for  the  latter  case  are  easily  obtained  from  those  for  the  former 
case  in  a  formal  way  by  simple  and  systematic  substitutions  that  have  been  given  explicitly. 

The  rapidity  with  which  a  particular  computational  program  converges  is  largely  determined 
by  the  abundance  in  the  mixture  of  those  constituents  chosen  as  components  relative  to  the  abun¬ 
dance  in  the  mixture  of  the  remaining  derived  constituents.  The  equilibrium  composition  of  car¬ 
bon,  hydrogen,  oxygen,  nitrogen  systems  varies  greatly  if  considerable  variation  of  the  gross  com¬ 
position  parameters  and  of  the  temperature  and  pressure  are  to  be  considered.  Consequendy,  the 
program  that  is  suitable  for  a  given  set  of  conditions  of  state  and  composition  may  converge  very 
slowly  or  even  diverge  if  applied  to  another  region  of  different  characteristics.  For  the  efficient 
consideration  of  all  possible  gross  compositions  over  considerable  ranges  of  temperature  and 
pressure  a  number  of  different  programs,  corresponding  to  different  choices  of  the  components, 
will  be  required.  Each  of  these  programs  will  be  most  appropriately  employed  for  particular  re¬ 
gions  of  the  complete  set  of  systems.  It  is  unnecessary  to  specify  closely  the  region  of  appli¬ 
cability  of  a  particular  program,  since  there  is  considerable  overlap  with  the  regions  appropriate 
for  alternative  programs.  It  will  be  noted  that  the  procedure  for  setting  up  a  particular  program  is 
routine,  once  the  constituents  to  be  considered  have  been  specified  and  a  choice  of  components 
has  been  made. 

It  will  generally  be  assumed  that  the  equilibrium  composition  of  the  combustion  products  of 
fuel-oxidant  systems,  composed  only  of  compounds  of  carbon,  hydrogen,  oxygen,  and  nitrogen, 
can  be  adequately  represented  by  a  consideration  of  the  constituents:  C( solid)’  CO,  C02r  ^2’ 
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H2O,  N2,  02,  O,  OH,  H,  NO,  N,  NH^,  CH4.  Solid  carbon  will  be  assumed  to  be  in  the  form  of 
graphite.  If  the  oxidant  is  air,  argon  will  be  introduced  as  an  additional  constituent.  Argon  will 
necessarily  be  chosen  as  an  additional  component,  and,  as  it  is  inert,  it  does  not  affect  in  any 
important  way  the  computational  programs.  We  assume  that  the  ideal  gas  law  provides  an  ade¬ 
quate  equation  of  state.  In  the  previous  section,  the  modifications  necessary  to  take  account  of 
gas  imperfections  have  been  described 

The  Generol  Method  for  Lean  Mixtures 


Fuel-oxidant  mixtures  containing  an  excess  of  oxygen  over  the  amount  required  for  stoichio¬ 
metric  conversion  of  the  carbon  to  carbon  dioxide  and  the  hydrogen  to  steam  are  called  lean  mix¬ 
tures  and  are  characterized  by  C>  1/2,  where  C  is  the  composition  parameter  defined  by  equation 
(3.3).  The  predominant  constituents  of  the  combustion  products  satisfying  the  condition  of  inde¬ 
pendence  required  for  dieir  choice  as  components  are 

C02,  up,  O*  N2, 


except  at  very  high  temperatures.  A  more  convenient  set  of  components  for  mixtures  at  such  tem¬ 
peratures,  (c  a.  4000°  K.)  is  described  in  a  later  section.  If  the  oxidant  is  air,  a  fifth  component, 
argon,  must  be  considered. 

The  solution  of  equations  (3.13)  and  (3.14)  for  this  choice  of  components  results  in  the 
equations 

_ 4  (C/O) _ 

q  (<“°2)  _  2  +.(H/0)  +  2  (N/O)  +  4  (A/O)  ’ 

_ 2(H/Q) _ 

q  (H20)  -  2  +  (H/G)  +  2  (N/o)  +  4  ( A/O)  ’ 

_  2  -  4  (C/O)  -  (H/O)  (5  l} 

q(°2)  “  2  +  (H/O)  +  2  (N/O)  +  4  (A/O)  ’ 

_  _ 2  (N/O) _ 

q  (n2)  -  2  +  (H/O)  +  2  (N/O)  +  4  (A/O)’ 

4  (A/O) 

q  (A)  _  2  +  (H/O)  +  2  (N/O)  +  4  (A/O)  ’ 


relating  the  normalized  stoichiometric  constants  to  the  atomic  ratios.  These  relations  can  be 
employed  together  w'ith  equations  (3.6)  or  (3-7)  to  evaluate  the  stoichiometric  constants  for  as¬ 
signed  values  of  the  parameters  descriptive  of  the  gross  composition  of  the  system.  If  the  oxi¬ 
dant  is  not  air,  (A/O)  =  0,  q  (A)  =0,  and  we  obtain 
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( 1-A)  (1-B)(1-Q 

q  ( C°2)  "( 1-C)  ( 2A  +  B  -  AB)  +  C  ( 1-A)  ’ 

2B  ( 1-A)  ( 1-C) 

q  (H2°)  ~(  i-C)  ( 2A  +  B  -  AB)  +  C  ( 1- A)  ’ 

_ (1-A)  (2C-1) _  (5'2) 

q  (°2)  ~  ( 1-C)  (  2A  +  B  -  AB)  +  C  ( 1-A)  ’ 

_ 2A  ( 1-C) _ 

q  (N2)  "(1-C)  (2A  +  B  -  AB)  +  C  ( 1-A)  ’ 

by  substitution  of  equations  (3.5)  into  equations  (5.1). 

In  lean  mixtures,  we  assume  that  methane  and  ammonia  are  absent  from  the  equilibrium  mix¬ 
ture,  an  assumption  that  is  confirmed  by  test  calculations  for  the  temperatures  of  interest  in  com¬ 
bustion  research.  The  chemical  reactions  in  the  form  of  equations  (3.  11)  for  the  formation  of  the 
dependent  constituents  from  the  components  are 

CO  2  -( 1/  2)0 2  =  CO, 

29  -( 1/ 2)0 2  =H2, 

(1/2)0  2  =0, 

(1/2)1120  +(1/4)0  2  =  01! 

(1/2)1120  -(1/4)0  2  =  «, 

( 1/  2)N  2  +(1/2)0  2  =NO, 

( 1/  2)N  2  =N, 

in  which  the  symbols  of  all  of  the  components  are  written  on  the  left-hand  side  of  the  equations. 
The  coefficients  Vj:  of  these  equations  can  be  summarized  in  matrix  form, 

j  C02  H2O  02  N2  A 

0-1/2  0  0 
1  -1/2  0  0 

0  1/2  0  0  (5.3) 

1/2  1/4  0  0 

1/2  -1/4  0  0 

0  0  1/2  0 

0  0  1/2  0 


(vij)  =  CO  1 

H  2  0 

O  0 

OH  0 

H  0 

NO  0 

N  0 
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The  constants  kj  of  the  mass  action  expressions,  equations  (4.9),  are  given  for  stated  values  of 
the  temperature  and  pressure  by 


k(CO)  =  p-  1/2K(CO), 

k(H)  =  p*3/4K(H), 

k(H2)  =p-1/2K(H2), 

k(NO)  =  K(NO), 

(5.4) 

k(O)  =p-l/2K(0)) 

k(N)  =  p"l/2K(N), 

k(OH)  =p‘  !/4k(OH), 

where  the  thermodynamic  equilibrium  constants 

can  be  calculated  from  the  equilibrium 

constants 

of  formation  by  means  of  the  relations 

K(CQ)  =  Kf(CO)/Kf(CD2), 

K(H)  =  Kf(H)/  [Kf(H20)]  V2, 

K(H2)  =  l/Kf(H20), 

X(NO)  =  Kf(NO), 

(5.5) 

K(O)  =  Kf(O), 

K(N)  =  Kf(N), 

K(OH)  =  Kf(OH)/  [Kf(H20)]  l/2. 


In  writing  down  equations  (5.5),  we  have  taken  account  of  the  fact  that  the  equilibrium  constant 
of  formation  of  an  element  is  unity. 

The  mole  fractions  of  the  dependent  constituents  are  determined,  for  assumed  values  of  the 
mole  fractions  of  the  components,  by  equations  (4.9)  in  the  form 


x(CO)  =  k(C0)x(C02)/  [x(02)]  1/2, 
x(ll2)  =  k(Il2)x(H20)/  [x(02)]  l/2, 
x(0)  =  k(0)[x(02)]  V2f 

x(OH)  =  k(OH)  [x(H20)J  1'/2[x(02)]  !/4 ,  (5.6) 

x(H)  =  k(H)  [x(H20)j  V2/  [X(02)l  1/4, 
x(NO)  =  k(NO)  [x(02)]  1/ 2  [ x( N 2) ]  V-2, 
x(N)  =  k(N)  [x(N2)]  V2. 


The  values  of  the  mole  fractions  of  the  components  leading  to  an  exact  calculadon  of  the  equi¬ 
librium  composition  satisfy  the  conditions 

F(C02)  =0,  F(ll20)  =0,  F(02)  =0,  F(N2)  =G, 

where  the  functions,  defined  by  equations  (4.11),  are  given  by 
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F(C02)  =  Gq( CO 2)  -  X(C02)  -  x(CO), 

FfHjO)  =  Gq(H20)  -  -  x(H2)  -  x(OH)/2-  x(H)/2, 

F(02)  =  Gq(02)  -  x(02)  +  x(CO)/2  +  x(H2)/2  -  x(0)/2  (5.7) 

x(OH)/4  +  x(H)/4  -  x(NO)/2 , 

F(N2)  =  Gq(N2)  -  x(N2)  -  x(NO)/2-  x(N)/2, 

and  where 

G  =  1  -  x(CO)/2  -  x(H2)/2  -  x(0)/2  -  x(OH)/4  -  3x(H)/4  -  x(N)/2.  (5.8) 


The  application  of  the  Newton-Raphson  method  to  the  improvement  of  an  approximate  set  of 
values  of  the  mole  fractions  of  the  components  requires  the  solution  of  equations  (4.  12)  in  the 
form 


r"'.  .  .  .  N  _ >.  __  -S 


All  a12  a13  a14 

h(C02) 

= 

F(G02) 

a21  a22  a23  A24 

hCHp) 

F(H20) 

a31  a32  a33  a34 

h(02) 

f(o2) 

^41  a42  a43  A44. 

h(N2) 

F(N2) 

where  h(C02t,  h(H20),  h(02),  and  h(N2)  are  estimates  of  the  fractional  errors  of  a  particular 
approximation  to  the  mole  fractions  of  the  respective  components,  to  be  employed  as  indicated 
by  equations  (4.  13)  in  obtaining  the  next  approximation,  and  where  the  coefficients,  defined  by 
equations  (4. 14),  are  given  by 


An  =  un-  Viq(C02), 
A12  =  Uf2  -  V2q(C02), 
A13  =  u13  '  v3<l(G02), 
A14  =  ui4  -  V4^C02), 
a21  =  u12  *  v1^h2°)» 
a22  =  u22  *  v2q(H20), 
A23  =  ^23-  V3q(H20), 
a24  =  u24  '  v4^h2P)> 


a31  =  u13  •  viq(°2). 
a32  =  u23  '  v2q(02)> 
A33  =  U33  *  v3q(°2). 

A34  =  U34  -  .V4q(02), 
A41  =  U14-  V^Nz), 
a42=  U24- V2q(N2), 
A43  =  U34  -  V3q(N2), 
A44  =  U44  -  V4q(N2), 
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wi  th 


Un  =  x(C02)  +  x(CO),  U12  =  O,  U13  =  -  x(CO)/2, 
u14  =  O,  U22  =  x(H20)  +  x(H2)  +  x(OH)/4  +  x(H)/4, 

U23  =  x(H2)/2  +  x(OH)/8-  x(H)/8,  7  U24  =  O, 

U33  =  x(02)  +x(CO)/4  +x(H2)/4  +x(0)/4  +x(OH)/16  +x(H)/16  +  x(NO)/4, 
U34  =  x(NO)/4,  U44  =  x(N2)  +x(NO)/4  +  x(N)/ 4, 


Vx  =  -x(CO)/2,  V2  =  -x(H2)/2-  x(OH)/8-  3x(H)/8, 

V3  =  x(CO)/4  +  x(H2)/4-  x(0)/4  -  x(OH)/16  +3x(H)/l6,  V4  = -x(N)/4. 


If  the  oxidant  is  air,  the  value  of  the  mole  fraction  of  argon  in  the  equilibrium  mixture  can  be 
calculated  from 

x(A)  =  Gq(A).  (5.10) 

Ihe  total  number  of  moles  n  contained  in  MQ  grams  of  mixture  is  equal  to  G' where 

Mq  =44.010  q(C02)  +  18.016  q(H20)  +  32.000  q(C>2)  +  28.016  q(N2)  +  39-9^4  q(A), 

and  where  the  quantity  G  is  defined  by  equation  (5.8). 

A  Special  Method  for  Lean  Mixtu  re  s 


The  special  method  of  computation  can  profitably  be  employed  in  calculations  for  conditions 
such  that  the  convergence  of  the  general  method  is  inconveniently  slow.  This  will  be  the  case 
for  mixtures  near  the  stoichiometric  point  C  =  1/ 2  and  for  mixtures  at  temperatures  such  that  dis¬ 
sociation  of  the  components  into  atomic  species  is  important.  Although  the  special  method  con¬ 
verges  more  rapidly  than  the  general  method,  each  iteration  of  the  special  method  involves  a  con¬ 
siderably  longer  and  less  systematic  computational  routine  than  is  required  by  die  general 
method. 

Vfe  take  the  variable  (_l,  defined  by  equation  (4.17),  to  be 

|i  =  [n(02)]  X/2  (5.11) 

and  employ  equation  (5.11)  to  eliminate  n(02)  from  the  set  of  equations  describing  the  equilib¬ 
rium  conditions.  The  stoichiometric  conditions,  equations  (4.1),  on  the  components  CO  j,  K20, 
and  N 2  are  given  by 
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n(C02)  +  n(CO)  =q(C02), 

n(H20)  +  n(H2)  +  n(OII)/2  +  n(H)/2  = 

n(N2)  +  n(NO)/2  +  n(N)/2  =q(N2). 

The  mass  action  expressions,  equations  (4.  19),  can  be  written  explicitly 


where 


n(CO)  =  k  '(CO)x/n, 

n(H2)  =  k'(H2)y2/h, 

n(O)  =k'(0)H, 
n(OH)  =  k'(OH)yM.I/2, 

k'(CO)  =  n 2  k( CO) , 
k'(li2)  =  n1/2k(ll2), 
k'(O)  =n1/2k(0)) 

k  '(OH)  =  n1/4k(OH), 


n(H)  =k'(H)y /u1/2, 
n(NO)  =  k  '(NO)z|i, 
n(N)  =  k  '(N)z, 

k  '(H)  =n3/4k(H), 
k  '(NO)  =  k(N'O), 
k  '(N)  =  n1/2  k(N), 


and  where  we  have  employed  the  abbreviations 

x=n(C02),  y2=n(H20),  z2  =n(H2). 

Substituting  equations  (5.13)  into  equations  (5. 14),  we  obtain 

x[l  +  k'(CO)/n]  -  q(C02)  =  O, 

y 2 [  1  +  k'(H2)/(i]  +  yiVCOH)!!1/2  +  k'(H)/li1/2]/2.  q(H20)  =  O, 
z2  +  z  [  k'(NO)|l  +  k'(N) ]/ 2  -  q(N2)  =  O, 


with  the  solutions, 


q(co2) 

1+  k'«D)/u 


y 


k  '(OH)  [1 1/7  2  +k'(H)/|i1/2 
4  +  4  k  '(H2)/ |i 


K'(OH)iil/2  +  k'(H)li1/2 


L  4  +  4  k  '(H2)/ |i  J 


+  q(H2o) 
l+k'(H2)/u’ 


Z 


k'(NO)|i  +k'(N)  + 
4 


'(NO)|i  +k'(N) 
4 


2 

+  q(N2). 


(5.12) 


(5. 13) 


( 5. 14) 


(5.15) 


(5.16) 
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Equations  (5.16),  together  with  equations  (5.11),  (5.13),  and  (5.15),  determine  the  equilibrium 
composition  for  assumed  values  of  n  and  [!. 

The  values  of  n  and  |1  leading  to  an  exact  calculation  of  the  equilibrium  composition  satisfy 
the  conditions, 


F(b')  =0,  F(n)  =0, 

where  the  functions  are  defined  by  equations  (4.22)  and  have  the  explicit  form 

F(U)  =q(Oo)  -  n(()?)  +  n(CO)/ 2+n( Hi)/ 2-n(0)/ 2-n(OH)/4+n(H)/4-n(NO)/ 2, 

(5.17) 

F(n)=n-l-n(C0)/2-  n(H2)/2-  n(0)/2-  n(OH)/4  -  3n(li)/4  -  n(N)/2. 

If  F  and  G  do  not  vanish,  the  application  of  the  Newton-Raphson  method  to  obtain  an  im¬ 
proved  pair  of  values  of  n  and  (i  requires  the  solution  of  equations  (4.23), 


"a,, 

a12 

h(u) 

= 

F(u) 

^A21 

a22. 

h(n) 

F(n), 

where  h(p)  and  h(n)  are  estimates  of  the  fractional  error  of  the  approximation  to  pi  and  n,  respec¬ 
tively,  to  be  employed  as  indicated  by  equations  (4.24)  in  obtaining  the  next  approximation,  and 
where  the  coefficients,  defined  by  equations  (4.25),  are  given  by 

An  =  2n(02)  +  n(C0)/2  +n(H2)/2  +  n(D)/2  +  n(OH)/8  +  n(H)/8  +  n(NO)/2 

+  U  [-n(C0)/2]  +V  [-  n(H2)  +  n(OH)/4  -  n(H)/4  ]  +  W  [n(NO)/2  ] 

A12  =  -n(C0)/4  -  n(H2)/4  +  n(0)/4  +  n(OH)/  16  -  3n(H)/  16 

+  X  [-  n(C0)/2]  +  Y  [-  n(H2)  +  n(OH)/4  -  n(H)/4]  +  Z  [n(N0)/2] 

A21  =  -  n(CO)/2-  n(H2)/2  +n(0)/2  +n(OH)/8-  3n(H)/8 

+  U  [n(CO)/2]  +v[n(H2)+n(OH)/4  +  3n(H)/4]  +W  [n(N)/2] 

A22  =  -  n  +  n(C0)/4  +  n(H2)/4  +  n(0)/4  +  n(OH)/ 16  +9n(H)/l6  +  n(N)/4 

+  X  [n(G0)/2]  +  Y  [n(fl2)  +  n(OH)/4  +  3n(Il)/4]  +  Z  [n(N)/2] 

where 

n(CO)  n(CO)/2 

U  _  q(C02)’  X  q(C02)  ’ 


n(H2)  -  n(0H)/4  +  n(H)/4 
q(fl20)  +  n(H2Q)  +  n(H2)  ’ 


n(H2)/2  +n(0H)/8  +3n(H)/8 
qtH^)  +n(H2G)  +  n(H2) 


n(  NO)/ 2 
q(N2)  +  n(N2)  ’ 


n(N)/4 

q(N2)  +  n(N2)  ' 
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If  the  oxidant  is  air,  the  concentration  of  argon  is  given  by 


n(  A)  =  q(  A).  (5.19) 

The  equilibrium  composition  is  expressed  in  terms  of  mole  fractions  by  dividing  the  number  of 
moles  of  each  constituent  by  the  total  number  of  moles  n  of  the  mixture.  The  mass  of  the  mix¬ 
ture  in  grams  containing  a  total  number  of  gram  moles  n  is  MQ,  where 

Nlo  =44.010  q(C02)  +  18.016  q(H20)  +  32.000  q(02)  +  28.016  q(N2)  +  39.944  q(A). 

The  General  Method  for  Rich  Mixtures  Without  Solid  Carbon 


Fuel-oxidant  mixtures  containing  a  deficit  of  oxygen  compared  to  the  amount  re  quired  for 
stoichiometric  conversion  of  the  carbon  to  carbon  dioxide  and  the  hydrogen  to  steam  are  called 
rich  mixtures  and  are  characterized  by  C<l/2,  where  C  is  the  composition  parameter  defined  by 
equation  (3.2).  If  C>Cj,  where  Cj  is  the  critical  value  of  C  defining  the  solid  carbon  boundary 
surface,  equation  (3.31),  the  equilibrium  mixture  will  consist  of  a  gas  phase  only.  The  predom¬ 
inant  constituents  of  the  combustion  products  satisfying  the  conditions  of  independence  required 
for  their  choice  as  components  are 

CO,  H2,  Hp,  N 2, 


except  at  very  high  temperatures.  A  more  convenient  set  of  components  for  mixtures  at  such  tem¬ 
peratures  (ca.  4000°  K.)  is  described  in  a  later  section.  If  the  oxidant  is  air,  a  fifth  component, 
argon,  must  be  considered. 

The  solution  of  equations  (3.13)  and  (3.14)  for  this  choice  of  components  results  in  the 
equations 

_ 2(C/Q) _ 

q(C°)  -  2(cSo)  +  (h/O)  +  (N/O)  +  2(A/0)  ’ 


q(H2) 


_ (H/O)  +  2(C/Q)  -  2 _ 

2 (C/O)  +  (H/O)  +  (N/O)  +  2(  A/O)  ’ 


_ 2-  2(C/Q) _ 

q( h2°)  -  2(c/0)  +  (H/O)  +  (N/O)  +  2(A/0)  ’ 


(5.20) 


_ (N/O) _ __ 

q(N2)  ”  2 (C/O)  +  (H/O)  +  (N/O)  +  2( A/O)  ’ 


q(  A) 


_ 2(4/0) _ 

2(  C/O)  +  (H/O)  '+  (N/O)  +  2 (A/O)  ’ 
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relating  the  normalized  stoichiometric  constants  to  the  atomic  ratios.  These  relations  can  be 
employed,  together  with  equations  (3.6)  or  (3.7),  to  evaluate  the  stoichiometric  constants  for  as¬ 
signed  values  of  the  parameters  descriptive  of  the  gross  composition  of  the  system.  If  the  oxi¬ 
dant  is  not  air,  (A/O)  =0,  q(A)  =0,  and  we  obtain 


q(  CO) 


(l-AMl-B) 

1  +  A  +  B  -  AB  ’ 


-  (1'A)  (1-Q  -  2C] 

q(l!2)  ”  (1-C)  (1  +  A  +  B  -  AB)’ 


q(Hp) 


(l-A)  [2C-  (l-B)  (1-C)] 

(  l-C)  (1  +  A  +  B  -  AB)  ’ 


q(N2) 


2A 

1  +  A  +  B  -  AB 


(5-21) 


by  substitution  of  equations  (3.5)  into  equations  (5.20). 

The  chemical  reactions  in  the  form  of  equations  (3.  11)  for  the  formation  of  the  dependent  con' 
stituents  from  the  components  are 


CO  +  up  -  H2  =  co2, 
21120  -  211 2  =  O  2 , 


h2o  -  H2  =  O, 

H20  -  Ih2  =  OH, 

2 


2 


H2O  +  I  N2  -  H2  =  NO  , 


CO  +  3H 2  -  H20  =  CH4, 
1  H2  +  I  N2  =  nh3, 


in  which  the  symbols  of  all  of  the  components  are  written  on  the  left-hand  side  of  the  equations. 
Ihe  coefficients  V q j  of  these  equations  can  be  summarized  in  matrix  form, 
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The  mole  fractions  of  the  dependent  constituents  are  determined,  for  assumed  values  of  the 
mole  fractions  of  the  components,  by  equations  (4.9)  in  the  form 


x(co2)  = 

k(C02)  x(C0)x(H20)/x(H2), 

x(02)  = 

k(02)  [ x ( H 20) ]  2/[x(H2)]  2, 

x(O)  = 

k(0)  x(Hp)/x(H2), 

x(OH)  = 

k(OH)x(H20)/[x(H2)]  1/2, 

x(H)  = 

k(H)  [ x ( H 2) ]  1/2; 

x(NO)  = 

k(NO)x(H20)[x(N2)]  1/2/x(H2)  , 

x(N)  = 

k(N)[x(N2)]  ^2, 

x(CH4)  = 

k  (CH4)x(CO)  [ x ( 1 1 2) ]  3/x(H2p), 

x(NH3)  = 

k(NH3)[x(H2)]  3/2[x(N2)]  l/2. 

The  values  of  the  mole  fractions  of  the  components  leading  to  an  exact  calculation  of  the  equi¬ 
librium  composition  satisfy  the  conditions 

F(CO)  =0,  F(H2)  =0,  F(H20)  =0,  F(N2)  =0,  (5.26) 


where  the  functions,  defined  by  e-quations  (4.  11),  are  given  by 

F( GO)  =  Gq ( CO)  -  x(CO)  -  x(C02)  -  x(CH4>, 

F(H2)  =  Gq(H2)  -  x(H2)  +  x(G02)  +  2x(02)  +  x(0)  +  x(OH)/2 
-x(H)/2  +  x(NO)  -3x(CH4)  -  3x(NH3)/2, 

(5.  27) 

F(H^O)  =  Gq(H20)  -  xCHjO)  -  x( Cj02)  -  2x(02)  -  x(0)  -  x(OH) 

-x(NO)  +  x(CH4), 

F(N2)  =  Gq(N2)  -  x(N2),-  x(NO)/2  -  x(N)/2  -  x(NH3)/2, 


and  where 


G  =  1  -  x(02)  -  x(0)  -  x(OH)/2-  x(H)/2  -  x(NO)/2  -  x(N)/2 
+  2x(CH4)  +  x(NH3). 


(5.28) 


The  application  of  the  Newton-Raphson  method  to  the  improvement  of  an  approximate  set  of 
values  of  the  mole  fractions  of  the  components  requires  the  solution  of  equations  (4.  12)  in  the 
form 
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s' 

A11 

a12 

a13 

a14 

"h(CO)V 

= 

F(GO) 

A2 1 

a22 

a23 

a24 

h(H2) 

F(H2) 

a31 

A32 

a33 

a34 

h(Hp) 

F(H20) 

a41 

\ 

a42 

a43 

a44 

h  (N2) 

F(N2) 

(5.29) 


where  h(CO),  h(H2),  h(H20),  and  h(N2)  are  estimates  of  the  fractional  errors  of  a  particular  ap¬ 
proximation  to  the  mole  fractions  of  the  respective  components,  to  be  employed  as  indicated  by 
equations  (4.  13)  in  obtaining  the  next  approximation,  and  where  the  coefficients,  defined  by  equa¬ 
tions  (4.  14),  are  given  by 


An  =  Un  •  viq(c°)> 

A12  =  U12  '  V2q(CO), 
A13  =  u  13  *  v3<l(CO), 
a14  =  L!14  -  V4q(CO), 

a21  =  u12  '  viq(H2)> 
a22  =  u22  ‘  v2q(H2>. 

a23  =  u23  '  v3q(H2>- 
a24  =  u24  *  V4q(H2), 

with 


a31  =  U13  -  Viq(Hp), 
a32  =U23-  V2q(H2p), 
A33  =  u33  ‘  VjqCHjO), 
A34  =  1-34  -  V4q(H2P), 

a41  =U14- Viq(N2), 
a42  =  ^24  *  V2q(N2), 
A43  =  U34  -  V3q(N2), 
A44  =  U44  -  V4q(N2), 


Uu  =  x(CO)  +  x(G02)  +  x  ( Cll  4) ,  U12  =  -x(C02)  +  3x(CH4), 

U i3  =  x ( CO 2)  -  x(Cll4),  L'14  =  O, 

U22  =  x(H2)  +  x(C02)  +  4x(02)+  x(O)  +  x(OH)/4  +  x(H)/4  +  x(NO)  +  9.x(CH4) 
+  9x(NH3)/4, 

U23  =  x(C02)  -  4x(02)  -  x(O)  -  x(OH)/2  -  x(NO)  -  3x(CH4), 

U24  =  -x(NO)/ 2  -  3x(NH3)/4, 

U33  =  x(H20)  +  x(C02)  +  4x(02)  +  x(O)  +  x(OH)  +  x(NO)  +  x(CH4), 

U34  =  x(NO)/2,  U44  =  x(N2)  +  x(NO)/4  +  x(N)/4  +  x(NH3)/4, 
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and 


Vx  =  2x(CH4), 

V2  =  2x(02)  +x(0)  +x(0H)/4  +x(H)/4  +  x(NO)/2  +6x(CH4)  +  3x(NH3)/2, 

V3  =  -2x(02)  -  x(O)  -  x(0H)/2-  x(NO)/2-  2x(CH4), 

V4  =  -x(NO)/4  -  x(N)/4  +  x(NH3)/4. 

If  the  oxidant  is  air,  so  that  q(A)  7  O,  the  value  of  the  mole  fraction  of  argon  in  the  equilib¬ 
rium  mixture,  can  be  calculated  from 


x(A)  =Gq(  A).  (5.30) 

The  total  number  of  moles  n  contained  in  MQ  grams  of  mixture  is  equal  to  G‘  where 

M  =28.010  q(CO)  +  2.016  q(II2)  +  18.016  q(H20)  +  28.016  q(N2)  +  39.944  q(A), 

and  where  the  quantity  G  is  defined  by  equation  (5.28). 

The  rate  of  convergence  of  this  method  of  computation  can  be  increased  substantially  if  the 
mole  fractions  of  dependent  constituents  are  calculated  from  an  approximate  set  of  values  of  the 
variables,  y^,  y2,  y3,  and  x(N2),  where 


y  i  -x(co)  t  x(co2), 

y2  =x(H2)  -  x(C02),  (5.31) 

y3  =x(h2o)  +  x(co2). 


We  may  determine  a  value  of  x(C02),  consistent  with  given  values  of  yj,  y2,  and  y3  and  with  the 
first  of  equations  (5.26),  by  means  of 


c(C02)  =  -b  +  \|b2  +  c  for  k(C02)  <  1, 
- - 

c ( CO 2)  =  -b  -  \  b2  +  c  for  k(C02)  >  1, 


where 

.  Y2  +  (y  1  +  y3>  k(C02) 

2  ■  [  1  -  k(C02)]  '* 


where  k(C02)  -  1,  the  solution  is 


y  1  y3  k(C02) 
1  -  k  (C02) 


x(G02) 


y\V3 

yi  +  y2  +  y  3 


(5.32) 


(5.33) 
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The  values  of  x(CO),  x(H2),  and  x(H20)  can  then  be  determined  for  given  y  j,  y2,  and  y^  by  equa¬ 
tions  (5.31),  and  they  may  be  employed  with  equations  (5.26)  to  complete  the  determination  of  the 
equilibrium  composition  for  assumed  values  of  y  j,  y2,  yy  If  this  determination  does  not 

satisfy  the  criteria  for  convergence,  the  next  approximation  can  be  obtained  by  means  of 

y'l=  71  [1  +  h  (CO)]  +  x(C02)  [-h (H2)  +  h(H20)], 

y2  ~  y2  [1  +h(li2)]  +  x(C02)  [-h(CO)  +  2h(H2)  -  h(H20)]  ,  (5.34) 

y  3  =  y 3  [l  +  h  (HjO)]  +  x(C02)  [h(CO)  -  h(ll2)]  , 
x'(N2)  =  x(N 2)  [l  +  h (N2)]  , 


!  I  1 

where  h(CO),  h(H ^>),  h(H20),  h(N2)  are  the  solutions  of  equations  (S-29)  and  where  yj,  y2,  y3, 
x'(N2)  are  the  new  approximations  to  yj,  y2,  y3,  and  x(N2),  respectively. 


A  Special  Method  for  Rich  Mixtures  without  Solid  Carbon 

The  special  method  of  computation  can  profitably  be  employed  in  calculations  for  conditions 
such  that  the  convergence  of  the  general  method  is  inconveniently  slow.  This  will  be  the  case 
for  mixtures  near  the  stoichiometric  point  C  =  1/2  and  for  mixtures  at  temperatures  such  that  dis¬ 
sociation  of  the  components  into  atomic  species  is  important.  Although  the  special  method  con¬ 
verges  more  rapidly  than  the  general  method,  each  iteration  of  the  special  method  involves  a  con¬ 
siderably  longer  and  less  systematic  computational  routine  than  is  required  by  the  general 
method. 

At  temperatures  sufficiently  high  so  that  methane  and  ammonia  can  be  assumed  to  be  absent 
from  the  equilibrium  mixture,  we  take  the  variable  |i,  defined  by  equation  (4.17),  to  be 


M.  =n(H20)/n(Ii2V  (5.35) 

and  employ  equation  (5.35)  to  eliminate  n(ll20)  from  the  set  of  equations  describing  the  equilib¬ 
rium  conditions.  The  stoichiometric  conditions,  equations  (4. 1),  can  be  written  in  the  form, 


n(CO)  +  n(C02)  =  q(CO), 

n ( H 2)  [  1  +  |i]  +  n(0H)/2  +  n(H)/2  =  q'(H2)  =  q(H2)  +  q(H20),  (5.36) 

n(N  2)  +  n(NO)/2  +  n(N)/2  =  q(N2), 

where  the  second  equation  is  a  linear  combination  of  the  stoichiometric  conditions  on  the  com¬ 
ponents  11 2  and  H2O.  The  mass  action  expressions,  equations  (4.19),  can  be  written  explicitly 


n(02)  =  k'(02)  d2, 
n  (O)  =  k'(O)  d, 

n(NO)  =  k  '(NO)  z  d, 
n(N)  =  k'(N)z, 
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n(OH)  =k'(OH)yd, 
n  (H)  k  '(H)  y, 
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(5.37) 


where 


k  '( CD  2)  =  k(C02),  k'(02)  =  nk(02), 

k  '(OH)  =n1/2k(Ok!),  k'(O)  =nk(0), 

k'(H)  =  n 2k  (H),  k'(NO)  =  n1/2k(N0)> 

k'(N)  =  n1/2k(N), 

and  where  we  have  employed  the  abbreviations 

x  =  n(CO),  y2  =  n(H2),  z2  =  n(N2). 
Substituting  equations  (5. 37)  into  equations  (5.36),  we  obtain 

x  [l  +  k'(C02)|i.]  -q(CO)  =  O, 

y2(l  +  H)  +  y  [k'(OH)n.  +  k'(H)]/2-  q'(ll2)  =  0, 

z2  +  z  [k  '(NO)  |i  +  k'(N)]/2-  q(N2)  =  0, 

with  the  solutions 

-  q(°°) 

X  ~  1  +  k  '(C02)(_i  ’ 

k  '(OH)  H  +  k  '(H) 
y  '  4(1+  n) 


2 

+  q  ( N  2) . 

Equations  (5.40),  together  with  equations  (5.35),  (5.37),  and  (5.39),  determine  the  equilibrium 
composition  for  assumed  values  of  n  and  |_t. 

The  values  of  n  and  (i  leading  to  an  exact  calculation  of  the  equilibrium  composition  satisfy 
the  conditions, 

F(  |_l)  =0,  F(n)  =0, 

where  the  functions  are  defined  by  equations  (4.22)  and  have  the  explicit  form, 

F(ji)  =  q(H?0)  -  n(HTO)  -  n(C02)  -  2n(0?)  -  n(O)  -  n(OH)  -  n(NO), 

(5.41) 

F(n)  =  n  -  1-  n(02)  -  n(O)  -  n(OH)/2  -  n(H)/2-  n(NO)/2-  n(N)/2. 


k  '(NO) |i  +  k'(N) 
z  -  - -  + 


k  '(NO)  |i  +  k'(N)~ 


k  '(OH)H  +  k  '(H) 


4  ( 1  +  H)  J 


q'(H2) 

1  +  Lt  ’ 


(5.40) 


( 5.  38) 


(5.39) 
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If  F  and  G  do  not  vanish,  the  application  of  the  Newton-Raphson  method  to  obtain  an  im* 
proved  pair  of  values  of  n  and  (1  requites  the  solution  of  equations  (4.23), 


"a„ 

\ 
A 12 

s  "1 
h(n) 

= 

F(H) 

A21 

a22 

y 

h(n) 

F(n) 

<  > 

(5.42) 


where  h(|i)  and  h(n)  are  estimates  of  the  fractional  error  of  the  approximation  to  |i  and  n,  respec¬ 
tively,  to  be  employed  as  indicated  by  equations  (4.24)  in  obtaining  the  next  approximation,  and 
where  the  coefficients,  defined  by  equations  (4. 25),  are  given  by 

Aj  i  =  n  (H2O)  +  n(C02)  +  4n(C>2)  +  n(0)  +  n(OH)  +  n(NO) 

+  U  n(G02)  +  V  [2n(H20)  +  n(OH)]  +  W  n(NO), 

A12  =  2n  (02)  +n(0)  +n(OH)/2  +n(NO)/2 

+  Xn(G02)+Y  [2n(Hp)  +  n(OH)]  +  Zn(NO), 

A21  =  2n(02)  +n(0)  +n(0H)/2  +n(NO)/2 

+  V  [n (0H)/2  +n(H)/2]  +  W  [n(NO)/2  +n(N)/2], 

A22  =  -n  +n(02)  +n(0)  +n(OH)/4  +  n(H)  +n(NO)/4  +  n(N)/4 
+  Y  [n(0H)/2  +  n(H)/2]  +  W  [n(NO)/2  +n(N)/2], 


where 


n(G02) 

'q(co)  ’ 


n(H20)  -  n(OH)/2 
q'(H2)+  n(H2)  +n(H20)  ’ 


n  (OH)/ 4  +n(H)/4 
q'(H2)  +  n(H2)  +n(H20)  ’ 


n(NO)/2 
q(N2)  n(N2) 


n  (N0)/4  +n(H)/4 
q(N2)  +n(N2) 


If  the  oxidant  is  air,  the  concentration  of  argon  is  given  by 

n(A)=q(A).  (5.43) 

The  equilibrium  composition  is  expressed  in  terms  of  mole  fractions  by  dividing  the  number  of 
moles  of  each  constituent  by  the  total  number  of  moles  n  of  the  mixture.  The  mass  of  the  mix¬ 
ture,  in  grams,  containing  a  total  number  of  gram  moles  n  is  MQ,  where 

M  =28.010  q(CO)  +  2.016  q(H2)  +  18.016  q(H20)  +  28.016  q(N2)  +  39.944  q(A). 
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It  should  be  noted  that  this  computation  program  remains  determinate  for  the  case  where  the 
system  contain-s  no  hydrogen  (B  =  O),  since  |I  remains  finite  as  nfh^)  and  nfl^O)  tend  to  zero. 

A  Special  Method  for  Homogeneous  Mixtures  at  High  Temperatures 

At  high  temperatures  (ca.  4,000°  K.),  a  considerable  amount  of  dissociation  results  in  the  for¬ 
mation  of  atomic  species  at  the  expense  of  the  substances  employed  as  components  tor  the  pre¬ 
viously  described  computational  methods  for  both  lean  and  rich  mixtures.  Under  these  circum¬ 
stances,  these  methods  are  inconveniently  slow  in  converging,  and  it  is  useful  to  employ  a 
special  computation  program  for  these  temperatures.  The  predominant  constituents  of  the  com¬ 
bustion  products  satisfying  the  conditions  of  independence  required  for  their  choice  as  compo¬ 
nents  are 

CO,  11,  O,  N> 


If  the  oxidant  is  air,  a  fifth  component,  argon,  must  be  considered.  At  these  temperatures,  am¬ 
monia  and  methane  need  not  be  considered  as  possible  constituents  of  the  equilibrium  mixture. 


The  solution  of  equations  (3.13)  and  (3-14)  for  this  choice  of  components  results  in  the 
equations 

2(  C/O) 

q(CO)  = 


q(H)  = 

q(0)  = 
q(N2)‘= 
q(  A)  = 


2  +2(11/0)  +(N/0)  +  2( A/O)  ’ 
2(11/0) 

2  +2(11/0)  +(N/0)  +  2(A/0)  ’ 
2-  2  ( C/O) 

2  +  2(11/0)  +(N/0)  +  2  (A/O)’ 
(N/O) 

2  +  2(11/0)  +(N/0)  +  2(  A/O)  ’ 
2 (A/O) 

2  +  2(11/0)  +(N/0).+  2(  A/O)’ 


(5.44) 


relating  the  normalized  stoichiometric  constants  to  the  atomic  ratios.  These  relations  can  be 
employed,  together  with  equations  (3.6)  or  (3.7),  to  evaluate  the  stoichiometric  constants  for  as¬ 
signed  values  of  the  parameters  descriptive  of  the  gross  composition  of  the  system.  If  the  oxi¬ 
dant  is  not  air,  ( A/O)  =  0,  q(A)  =0,  and  we  obtain  by  substitution  of  equations  ( 3.  5) 


( 1-A)  ( 1-B)  ( 1-Q 


q(  co)  = 

q(H)  = 

q(0)  = 

q(N2)  = 


2C(  1-  A)  +  2(  1-C)  (A  +  2B  -  2AB)  ’ 
2B  ( 1-A  ( 1-Q 

C(l-A)  +  (  1-C)  (A  +  2B  -  2AB)  ’ 

(1-A)  (B  +  3C-  BC-  1) 
2C(1-A)  +  2(1-C)  (A  +  2B-  2AB)  ’ 
A  (1-C) 

C(l-A)  +  ( 1-C)  (A  +  2B  -  2AB)  ' 


(5.45) 
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The  stoichiometric  constants  for  this  program  can  be  computed  from  those  for  the  program  for  rich 
mixtures  without  solid  carbon  at  lower  temperatures  by  means  of  the  relations 

q(  CO)  =  q(CO)(rV  [q(C0)(r)  +  2q(H2)(r)  +  3q(H20)(r)  +  q(N2/r)]  , 

q(H)  =  [2q(H2)(r)  +  2q(Hp)(r)]/  [q(CO)(r)+2q(H2)(r)  +  3q(H20)(r)  +  q(N2)(r)], 

(5.46) 

q(0)  =  q(Hp)(rV  [q(CO)(r)+2q(H2)(r)  +  3q(H20)(r)  +q(N2)(r>]  , 

q(N2)  =  q(N2)(rV[q(C0)(r)  +  2q(H2)(r)+3q(H20)(r)+q(N2)(r)], 

where  the  superscript  r  is  employed  to  designate  the  previously  defined  constants  for  homoge¬ 
neous,  rich  mixtures  at  lower  temperatures.  Similarly,  in  terms  of  the  constants  for  lean  mixtures 
at  lower  temperatures, 

q( CO)  =  q(C02)(^V  [2q(C02)<^)  +  3q(H20 +  2q(02)(^)  +  q(N2)(^)], 

q(H)  =  2q(H20)('ZV  [2q(C02)('/)  +  3q(H20)^>  +  2q(02)(^>  +  q(N2)(^)],  *  ' 

q(O)  =  [q(C02)(^)  +  q(Hp)(^)  +  2q(02)(^>]  /  [  2q(C02)(/>  +  3q(H20)(^)  +  2q(02)<^ 

+  q(N2)(^)], 

q(N2)  =  q(N2)(^)/  [2q(O02)('i’)  +  3q(H20)^  +  2q(02)(^>  +  q(N2)(^)]  , 

where  the  superscripts^  is  employed  to  designate  the  previously  defined  constants  for  lean  mix¬ 
tures  at  lower  temperatures. 

CO  +  o  =co2, 

2H  =H2, 

2H  +  O  =H20, 

2  0  =02, 

0  +  H  =OH, 

( 1/  2)N  2  +  O  =NO, 

( 1/  2)N  2  =N. 

The  coefficients  V  ••  of  these  equations  can  be  summarized  in  matrix  form, 
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(Vij) 


=  CO- 


H- 


Hp 

o2 

OH 

NO 

N 


CO 


H 


O 


1 

0 

0 

0 

0 

0 

0 


0 

2 

2 

0 

1 

0 

0 


1 

0 

1 

2 

1 

1 

0 


0 

0 

0 

0 

0 

1/2 

1/2 


(5.48) 


The  constants  Iq  of  the  mass  action  expressions,  equations  (4.9),  are  given  for  stated  values  of 
the  temperature  and  pressure  by 


k(C02)  =  pK(C02), 
k(H2)  =  pK(H2), 

k(H20)  =  p2K(Hp), 
k(02)  =  pIC(02), 


k  (OH)  =  pK(OH), 
k  (NO)  =  p^KfNO), 
k(N)  =p-l/2K(N), 


(5.49) 


where  the  thermodynamic  equilibrium  constants  can  be  calculated  from  the  equilibrium  constants 
of  formation  bv  means  of  the  relations 


K(C02)  =  Kf  (C02)/Kf  (CO)Kf(O), 

K(H2)  =l/[Kf(H)]2, 

K(H20)  =  Kf(H20)/[Kf(H)]  2  Kf(O), 

K(02)  =l/[Kf(0)]2. 

In  writing  equations  (5.50)  we  have  taken  account  of  the  fact  that  the  equilibrium  constant  of  for 
mation  of  an  element  is  unity. 

The  mass  action  expressions,  equations  (4.  19)  can  be  written  in  the  form 


n(C02)  =  k  '(C02)  xp, 

n  (OH)  =  k  '(OH)  yp, 

n(H2)  =k'(H2)y2, 

n  (NO)  =  k  '(NO)  zp, 

n(H20)  =  k'(H20)  y2p, 

n  (N)  =k'(N)z, 

n(02)  =  k'(02)p2. 

where  we  have  introduced  the  variable 


K(OH)  =  Kf(OH)/Kf  (H)Kf(O), 
K(NO)  =  Kf(NO)/Kf  (O), 

K(N)  =  Kf(N), 


and  the  abbreviations 


p  =n(0), 


(5.52) 


x  -n(CO),  y  -  n(H),  z2  -n(N2), 


(5.53) 
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and  where 


k'(C02)  =  n-1k(C02), 

k  '(OH)  =  n*  !k  (OH), 

k'(H2)  =n-lk(Il2), 

k  '(NO)  =  n'  !/2k  (NO), 

k'(li20)  =  n‘4(H20), 

k'(N)  =  n'1/2k(N), 

k'(02)  =  n“  !k(02). 

(5.54) 


We  employ  the  stoichiometric  conditions  in  the  form  given  by  equations  (4.  1).  The  conditions  on 
the  components  CO,  li,  and  N2  are 


n  ( CO)  +  n(C02)  =  q(CO), 

n(H)  +  2n(H2)  +  2n(H20)  +  n(OH)  =  q(H), 

n(N2)  +  n(NO)/2  +  n(N)/2  =  q(N2). 

Substituting  equations  (5.57)  and  (5.53)  into  equations  (5.55),  there  are  obtained 


(5.55) 


with  the  solutions 


x  [l  +  k'(CO)u]  -  q(CO)  -O, 

2y 2  [k'(H2)  +  k'(H20)|i]  +  y  [  1  +  k '(OH)|i]  -  q(H)  =0, 
z2  +  z  [k'(NO)(i  +  k'(N)]/2-  q(N 2)  =0, 


q(CO) 


1  +  k  '(CO)|l  ’ 


1  +  k'(OH)|i 

y  ’  *4k'(H2)  +  4k  '(H20)(i 


k  '(NO)[J.  +  k'(N) 


1  +  k'(OH)U 


4k  '(H2)  +  4k'(H20)|i 


q(H) 


2k '(H2)  +  2k'(H20).U  ’ 


k'(NO)p,  +  k'(N) 


q  (N2) . 


(5.56) 


Equations  (5.56),  together  with  equations  (5-51),  (5.52),  and  (5.53),  determine  the  equilibrium 
composition  for  assumed  values  of  n  and  |i. 

The  values  of  n  and  |J,  leading  to  an  exact  calculation  of  the  equilibrium  composition  satisfy 
the  conditions 

F(H)  =0,  F(n)  =0, 

where  the  functions,  defined  by  equations  (4.22),  are  given  by 

F(|i)  =  q(O)  -  n(O)  -  n(C02)  -  n(H20)  -  2n(02)  -  n(OH)  -  n(N0), 

.  (5.57) 

F(n)  =  n-1  +n(C02)  +n(H2)  +  2n(H20)  +  n(02)  +  n(OH)  +n(NO)/2-  n(N)/2. 
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If  F  and  G  do  not  vanish,  the  application  of  the  Newton-Raphson  method  to  obtain  an  im¬ 
proved  pair  of  values  of  n  and  p  requires  the  solution  of  equations  ( 4.  23), 


41 


l21 


42 


'22 


/  N 
h(u) 

= 

^F(u)^ 

h(n) 

4  / 

F(n) 

4  / 

(5.58) 


where  h(p)  and  h(n)  are  estimates  of  the  fractional  error  of  the  approximation  to  (1  and  n,  respec¬ 
tively,  to  be  employed  as  indicated  by  equations  (4.24)  in  obtaining  the  next  approximation,  and 
where  the  coefficients,  defined  by  equations  (4. 25),  are  given  by 


Au  =  n  (O)  +n(C02)  +n(H20)  +4n(02)  +n(OH)  +n(NO) 

+  Un(C02)  +  V[2n(H20)  +n(OH)]  +  Wn(NO), 

A12  =  -n(C02)  -  2n(H20)  -  2n(02)  -  n(OH)  -  n(N0)/2 
+  Xn(C02)  +  Y[2n(H20)  +  n(OH)]  +  Zn(NO), 

A21  =  - n ( CO 2)  -  2n(H20)  -  2n(02)  -  n(OH)  -  n(NO)/2 

+  U[-n(C02)]  +  v[-2n(H2)  -  4n(H20) -n  (OH)]  +  W[-n(N0)/2  +  n(N)/2], 
A22  =  -n  +  n(C02)  +  n(H2)  +  4n  (H20)  +  n(02)  +  n(OH)+  n  (N0)/4  +  n  (N)/4 

+  X[-n(C02)]  +  Y  [-2n(H2)  -  4n(H20)  -  n(OH)]  +  Z  [-n(N0)/2  +  n(N)/2]  , 


where 


n(C02) 
q(CO)  ’ 

2n(H20)  +  n  (OH) 
q(H)  +  2n(H2)  +  2n(H20)  ’ 

n(N0)/2 
q(N2)  +  n(N2) 


n  (C02) 
q(CO)  ’ 

n(H2)  +  2n(H20)  +.  n  (OH) 
q(H)  +  2n(H2)  +  2n(H20)  ’ 

n(NO)/2  -  n(N)/2 
q(N2)  +  n(N2) 


If  the  oxidant  is  air,  the  concentration  of  argon  is  given  by 


n(  A)  =q(A). 


(5.59) 


The  mass  of  the  mixture  containing  a  total  of  n  moles  is  MQ,  where 

M0  =28.010  q(CO)  +  1.008  q(H)  +  16.000  q(0)  +  28.016  q(N7)  +  39-944  q(A). 
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The  General  Method  for  Rich  Mixtures  with  Solid  Carbon 


Rich  fuel-oxidant  mixtures  that  are  characterized  by  C<C  p  where  C  is  the  composition  param¬ 
eter  defined  by  equation  (3.2)  and  is  the  critical  value  of  C  defining  the  solid  carbon  boundary 
surface,  equation  (3.31),  consist  of  a  gas  phase  in  equilibrium  with  solid  carbon.  The  predomi¬ 
nant  constituents  of  the  combustion  products  satisfying  the  conditions  of  independence  required 
for  their  choice  as  components  are 


* 


C(S)>  co>  h2>  n2* 

Because  the  system  is  heterogeneous,  it  possesses  only  two  degrees  of  freedom  for  the  specifi¬ 
cation  of  composition,  and  the  computation  program  is  based  upon  the  gaseous  components  only. 
If  the  oxidant  is  air,  a  fifth  component,  argon,  must  be  considered. 


The  solution  of  equations  (3.13)  and  (3.14)  for  this  choice  of  components  results  in  the 
equations 

2 

q(CO)  ~  2  +  (H/O)  +  (N/O)  +  2 (A/O)  ’ 


_ (H/O) _ 

q(H2)  -  2  +  (H/O)  +  (N/O)  +  2 (A/O)  ’ 

_ (N/O) _ 

q(N2>  “  2  +  (H/O)  +  (N/O)  +  2 (A/O)  ’ 


(5.60) 


q(  A)  = 


2(A/0) 


2  +  (H/O)  +  (N/O)  +  2 (A/O)  ’ 


relating  the  normalized  stoichiometric  constants  to  the  atomic  ratios.  These  relations  are  inde¬ 
pendent  of  the  amount  of  carbon  available,  which  is  assumed  to  be  in  excess.  They  can  be  em¬ 
ployed,  together  with  equations  (3.6)  or  (3.7)  to  evaluate  the  stoichiometric  constants  for  as¬ 
signed  values  of  the  parameters  descriptive  of  the  gross  composition  of  the  system.  If  the 
oxidant  is  not  air,  (A/O)  =0,  q(A)  =0,  and  we  obtain  by  substitution  of  equations  (3.5)  the 
relations 

q(CO)  =E(1-D), 

q(Il2)  =(1-D)(1-E),  (5.61) 

q(N2)  =d 

for  a  system  described  by  the  parameters  D  and  E. 

Vie  shall  determine,  in  addition  to  the  equilibrium  composition  of  the  gas  phase,  the  location 
of  the  solid  carbon  boundary  surface  in  the  form  =Bj(D,  E,  p,  T). 

The  chemical  reactions  in  the  form  of  equations  (3.11)  for  the  formation  of  the  dependent  con¬ 
stituents  from  the  components  are 
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2CO  -  C(s)  =  C02, 

CO  +  H2  -  C(S)  =  h2o, 
2CO  -  2<^s)  =  O, 

CO  -  c(s)  =  o, 

CO  +  j  H  2  -  C(  s)  =  OH, 

|h2  =  H, 

CO  +  ^  N  2  -  C(  s)  =  NO, 

jN2  =  N, 

9  s)  +  2H2=  OU- 

|h2  +  1n2  =  nh3. 


in  which  the  symbols  of  all  of  the  components  are  written  on  the  left-hand  side  of  the  equations. 
The  coefficients  v  jj  of  these  equations  can  be  summarized  in  matrix  form, 


C(s) 

CO 

h2 

n2 

A 

co2 

\ 

(-1) 

2 

0 

0 

0 

lip 

(-1) 

1 

1 

0 

0 

°2 

(-2) 

2 

0 

0 

0 

0 

(- 1) 

1 

0 

0 

0 

OH 

(- 1) 

1 

1/2 

0 

0 

H 

0 

0 

1/2 

0  . 

0 

NO 

(-D 

1 

0 

1/2 

0 

N 

0 

0 

0 

1/2 

0. 

ch4 

(1) 

0 

2 

0 

0 

NHr 

0 

0 

3/2 

1/2 

0 

D 

where  the  entries  in  the  column  for  solid  carbon  are  enclosed  in  parentheses  to  indicate  that 
these  elements  are  not  to  be  taken  into  account  in  summations  over  the  gaseous  components.  The 
constants  kj  of  the  mass  action  expressions,  equations  (4.9),  are  given  for  stated  values  of  the 
temperature  and  pressure  by 


* 
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k(C02)  =  pK(C02), 
k(H20)  =  pK(H20), 
k(02)  =  PK(02), 
k(O)  =  K(O), 
k  (OH)  =  p^^OH), 


k(H)  =p*1/2K(H), 
k  (NO)  =  p^KfNO), 
k(N)  =  p-V2x(N), 
k(CH4)=pK(CH4), 
k  (NH3)=pK(NH3), 


(5.63) 


where  the  thermodynamic  equilibrium  constants  can  be  calculated  from  the  equilibrium  constants 
of  formauon  by  means  of  the  relations 


K(C02)  =  Kf(C02)Kf(C)/  [Kf (CD)]  2, 
K(H20)  =  Kf(H20)Kf(C)/Kf(C0), 
K(02)  =  [Kf (C)]  2/  [Kf(CO>]  2 
K(O)  =  Kf  (0)Kf(C)/Kf(C0), 

K(Oll)  =  Kf(OH)Kf  (C)/Kf(CO), 


K(H)  =  Kf  (H), 

K(NO)  =  Kf  (NO)Kf(C)/Kf  (CO), 


K  (N)  =  Kf(N), 

K(CH4)  =  Kf  (CH4)/Kf(C), 
K(N113)  =  Kf(NH3). 


(5.64) 


In  writing  equations  (5.64),  we  have  taken  account  of  the  fact  that  the  equilibrium  constant  of  for¬ 
mation  of  a  gaseous  element  is  unity.  If  solid  carbon  is  assumed  to  be  in  the  form  of  graphite, 
then  Kf(C)  =  1.  However,  if  the  solid  carbon  is  assumed  to  exist  in  an  activated  form,  the  equi¬ 
librium  constant  of  formation  of  the  activated  form  can  be  used  explicitly  in  equations  (5.66). 

The  mole  fractions  of  the  dependent  constituents  are  determined  by  equations  (4.9)  in  the 
form 

x(C02)  =  k(CO)  [x(CO)]  2, 
x(H 20)  =  k(H20)  x(CO)  x(H2), 
x(02)  =  k(02)  [x(CO)]2, 

x(0)  =  k(O)  x(CO), 

x(OH)  =  k  (OH)  x(CO)  [x(ll2)]  1/2, 

r  m/7  (5'65) 

x(H)  =  k (H)  [x(H2)3  122 , 

x(NO)  =  k  (NO)  x(CD)  [x(N2)3  1/2> 

x(N)  =  k(N)  [x(N2)]  V2, 

x(ai4)  =k(cn4)  [x(h2)]  2, 

x(NH3)  =k(NH3)  [x(H2)]3/2  [x(N2)]  ^ 2 

for  assumed  values  of  the  mole  fractions  of  the  components.  The  values  of  the  mole  fractions  of 
the  gaseous  components  leading  to  an  exact  calculation  of  the  equilibrium  composition  of  the  gas¬ 
eous  phase  satisfy  the  conditions 


F(CO)  =0,  F(H2)  =Q,  F(N2)  =0, 
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where  the  functions,  defined  by  equations  (4.11)  are  given  by 


F(CO)  =  Gq(CO)  -x(GO)  -2x(C02)  -x(H20)  -2x(02)  -x(O)  -x(OH)  -x(NO), 

F(H2)  =  Gq(H2)  -x(H2)  -x(Hp)  -x(OH)/2-x(H)/2-2x(CH4)-3x(NH3)/2,  (5.66) 

F(N2)  =  Gq(N2) -x(N2) -x(NO)/2 -x(N)/2 -x(NH3)/2, 

and  where 

G  =  l+x(C02)  +  x(H20)  +  x(02)+x(0H)/2-x(H)/2+x(N0)/2-x(N)/2+x(ai4)+x(NH3).  (5.67) 

The  application  of  the  Newton-Raphson  method  to  the  improvement  of  an  approximate  set  of 
values  of  the  mole  fractions  of  the  gaseous  components  requires  the  solution  of  equations  (4.  12) 
in  the  form 


where  h(CO),  h(H2),  and  h(N2)  are  estimates  of  tiie  fractional  errors  of  a  particular  approximation 
to  the  mole  fractions  of  the  respective  gaseous  components,  to  be  employed  as  indicated  by  equa¬ 
tions  (4.  13)  in  obtaining  the  next  approximation,  and  where  the  coefficients,  defined  by  equations 
( 4. 14),  are  given  by 

AU  =  li  iq  -  Vjq(CO),  A23  =  U23  -  V3q(H2), 

A 12  =  U 12  *  V2  q(CO),  A3T  =  Uj3  -  V1q(N2), 

A13  =  u13  -  V3q(CO),  A32  =  u23  -  V 2  q ( N 2) ,  (5.69) 

a21  =  U21 '  vlcl(H2)>  a33  =  u33  *  v3<l(N2)> 

a22  =  u22  '  v2  cl(H2^> 

with 

un  =  x(CO)  +  4x(C02)  +  x(ll20)  +  4x(02)  +  x(O)  +  x(OH)  +  x(NO), 
u12  =  x(H20)  +  x(OH)/2,  u13  =x(NO)/2, 

(5.70) 

U 22  -  x(H2)  +  x(H20)  +  x((.'H)/4+x(H)/4  +  4x  (  CH4)/4+9x  (NH3)/4, 

U23  =  3x(NH3)/4,  U33  =x(N2)  +  x(NO)/4  +  x(N)/4  +  x(NH3)/4, 

and 

=  2x(C02)  +  x(H20)  +  2x(G2)  +  x(OH)/2  +  x(NO)/2, 

V2  =  x(II20)  +  x(OH)/4-  x(H)/4  +  2x(CH4)  +  3x(NH3)/2,  (5.71) 

V3  =  x(NO)/4  -  x(N)/4  +  x(NH3)/2. 
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If  the  oxidant  is  air,  so  that  q(A)  =0,  the  value  of  the  mole  fraction  of  argon  in  the  equilib¬ 
rium  mixture  can  be  calculated  from 

x(  A)  =  G  q(  A).  (5.72) 

The  total  number  of  moles  n  contained  in  MQ  grams  of  the  gaseous  mixture  is  G”^,  where 
Mo  =28.010  q(CO)  +  2.016  q(H2)  +  28.016  q(N2)  +  39.9  44  q(A) 

and  where  the  quantity  G  is  defined  by  equation  (5.67). 

In  order  to  determine  the  location  of  the  solid  carbon  boundary  surface,  the  number  of  gram 
atoms  of  carbon  contained  in  M  grams  of  gas  phase  in  the  form  of  compounds  of  carbon  is  cal¬ 
culated  by 

Q  =  —  [x(CO)  +  x(C02)  +  x(CH4)].  (5.73) 

C  G 


It  is  evident  from  equations  (5.66)  that  the  number  of  gram  atoms  of  hydrogen  in  the  same  weight 
of  gas  phase  is 


Qh  =  2q(H2). 


(5.74) 


If  these  relations  are  substituted  into  the  definition  of  B  given  in  equations  (3.3),  we  obtain 


q(H2) 

Bl  =  2Qc  +  q(H2)  • 


(5.75) 


Now  q( H 2)  is  a  function  of  D  and  E,  given  by  equations  (5.61),  and  Qq  is  a  function  of  the  equi¬ 
librium  composition,  and  therefore  an  implicitly  prescribed  function  of  D,  F,  the  temperature,  and 
the  pressure.  It  follows  that  equation  (5.75)  is  a  relation  of  the  form 


B  j  =B!(D,  E,  p,  T) 


for  the  location  of  the  solid  carbon  boundary  surface. 

If  the  amount  of  carbon  available  to  the  system  is  limited  and  specified  by  a  value  of  the  pa¬ 
rameter  B^Bj,  the  amount  of  solid  carbon  in  the  heterogeneous  mixture  at  equilibrium  is  given  by 
equation  (3.29),  which  can  be  written  in  the  form 

n(C(s))  =(1-D)  (1-E)(B1-B)/B1B,  (5.76) 

where  n(C^sp  is  the  number  of  gram  atoms  of  solid  carbon  in  equilibrium  with  MQ  grams  of  gas¬ 
eous  mixture. 

The  foregoing  computation  method  has  been  based  essentially  upon  a  specification  of  the  a- 
moun.'s  of  hydrogen,  oxygen,  and  nitrogen  available  to  a  system  with  an  excess  of  carbon.  This 
method  affords  the  most  convenient  basis  for  treating  systems  that  are  known  in  advance  to  be 
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heterogeneous.  However,  it  is  not  convenient  in  determining  the  range  of  applicability  of  a  com¬ 
putation  method  for  rich  homogeneous  mixtures  described  by  the  parameters  A,  B,  and  C.  For  the 
latter  purpose,  a  representation  of  the  solid  carbon  boundary  surface  in  the  form 

Ci  =C1(A,  B,  p,  T) 

is  desired,  requiring  a  computation  method  based  upon  a  specification  of  the  amounts  of  carbon, 
hydrogen,  and  nitrogen  available  to  the  system.  Such  a  description  of  the  composition  is  not  the 
natural  one  to  employ,  and  the  computation  program  for  this  case  can  not  be  formulated  by  the 
routine  application  of  the  general  principles.  However,  it  is  possible  to  devise  a  computation 
method  for  the  heterogeneous  mixture  whose  composition  is  described  by  specification  of  the  pa¬ 
rameters  A  and  B,  and  this  method  is  similar,  except  in  detail,  to  the  one  that  has  been  described. 

An  alternative  solution  of  equations  (3.  13)  and  (3.  14),  appropriate  to  the  present  application, 
can  be  written  in  the  form 

2(C/0) 

q(C0)  ”  2(C/0)  +  (H/O)  -  (N/O)  +  2 (A/O)  ’ 


_ (H/O) _ 

q(H2)  ”  2(C/0)  +  (H/O)  +  (N/O)  +  2(  A/ O)  ’ 

(N/O) 

q(N2)  ~  2 (C/O)  +  (H/O)  +  (N/O)  +  2 (A/O)  ’ 


(5.77) 


2(A/0) 

q(A)  ~  2 (C/O)  +  (H/O)  +  (N/O)  +  2( A/O)  ’ 


in  which  q(CO)  is  determined  by  the  relative  amount  of  carbon  available  instead  of  by  the  rela¬ 
tive  amount  of  oxygen  available,  as  in  equations  (5.60).  Equations  (5.77)  are  independent  of  the 
relative  amount  of  oxygen  available  to  the  system.  They  yield 


q(  CO) 


(l- A)  ( 1-B) 

1  +  A  +  B  -  AB  ’ 


2B(  1-  A) 

q(H?)  =  - ,  (5.78) 

M  1  1  +  A  +  B  -  AB  ’ 

W2A 

q(N2)  ”  1  +  A  +  B  -  AB  ’ 
for  a  system  described  by  the  parameters  A  and  B. 

The  mole  fractions  of  the  dependent  constituents  are  determined  by  equations  (5.65)  for  as¬ 
sumed  values  of  the  mole  fractions  of  the  components.  The  mass  action  constants  appearing  in 
these  expressions  are  defined  by  equations  (5.63)  and  (5.64). 
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The  values  of  the  mole  fractions  of  the  gaseous  components  leading  to  an  exact  calculation 
of  the  equilibrium  composition  of  the  gaseous  phase  satisfy  the  conditions 


F(CO)  =  O,  F(H2)  =  O,  F(N2)  =  O, 
where  the  functions  are  given  by 

F(CO)  =  Gq(CO)  -  x(CO)  -  x(C02)  -  x(CH4), 

F(H2)  =  Gq(H2)  -  x(H2)  -  x(H20)  -  x(OH)/2  -  x(H)/2  -  2x(CH4)  -3x(NH3)/2, 
F(N2)  =  Gq(N2)  -  x(N2)  -  x(NO)/2-  x(N)/2-  x(NH3)/2, 

and  where 


G  =  l-x(OH)/ 2  -x(H)/2  -x(NO)/2  -x(N)/2  +  2x(ai4)  +  x(NH3).  (5.80) 

The  first  member  of  equations  (5.79)  differs  from  the  corresponding  member  of  equations  (5.66) 
because  the  value  of  q(CO)  is  determined  for  the  former  case  by  the  relative  amount  of  carbon 
available  to  the  gaseous  phase  and  for  the  latter  case  by  the  relative  amount  of  oxygen  available 
to  this  phase.  As  a  result,  the  definition  of  G  given  by  equation  (5.80)  differs  from  that  given  by 
equation  (5.67). 

The  application  of  the  Newton-Kaphson  method  to  the  improvement  of  an  approximate  set  of 
values  of  the  mole  fractions  of  the  components  requires  the  solution  of  equations  (5.68),  the  coef¬ 
ficients  of  which  are  given  by  equations  (5.69).  However,  in  the  evaluation  of  equations  (5.69), 
the  definitions  of  the  quantities  Ujj  '  and  Vj,  given  for  the  former  case  by  equations  ( 5.70)  and 
(5.71),  respectively,  must  be  replaced  by  the  relations 

Un  =  x(CO)  +  x(COz)  +  x(CH4), 

U12=2x(CH4),  U13  =  O, 

U22  =  x(H2)  +  x(H20)  +  x(OH)/4  +  x(H)/4  +  4x(CH4)  +  9x(NH3)/4, 

U23  =  3x(NH3)/4,  U33  =  x(N2)  +  x(NO)/4  +  x(N)/4  +  x(NH3)/4 

and 


V  i  =  2x(CH4), 

v2  =  -x(OH)/4  -x(H)/4  +  4x(CH4)  +  3x(NH3)/2,  (5.82) 

V3  =-x(NO)/4  -x(.N)/ 4  +  x(NH 3)/ 2. 


If  the  oxidant  is  air,  the  value  of  the  mole  fraction  of  argon  is  given  by  equation  (5.72),  with 
G  defined  by  equation  (5.80).  The  total  number  of  moles  n  contained  in  MQ  grams  of  the  gaseous 
mixture  is  G*  where  WQ  is  again  given  by 

Mq  =28.010  q(CO)  +  2.016  q(H2)  +  28.016  q(N2)  +  39-944  q(A), 
and  where  the  quantity  G  is  defined  by  equation  (5.80). 
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The  number  of  gram  atoms  of  oxygen  contained  in  MQ  grams  of  gas  phase  is  given  by 

Q0  =  ^  tx(CO)  +  2x(C°2)  +  x(H20)  +  2x(02)  +  x(O)  +  x(OH)  +x(NO)]  .  (5.83) 

The  number  of  gram  atoms  of  carbon  and  hydrogen  in  the  same  weight  of  gas  phase  are 


Qc=q(CO),  Qpj  =2q(H2). 


(5.84) 


If  these  quantities  are  substituted  into  the  definition  of  C  given  in  equations  (3.2),  we  obtain 

i 


Cl  = 


2q(CO)  +  q(II2)  +  Qq" 


(5.85) 


Now'  q(CO)  and  q(H2)  are  functions  of  A  and  B,  given  by  equations  (5.78),  and  Qq  is  a  function 
of  the  equilibrium  composition,  and  therefore  an  implicitly  prescribed  function  of  A,  B,  the  tem¬ 
perature,  and  the  pressure.  It  follows  that  equation  (5.85)  is  a  relation  of  the  form 

C1  =  C|(A,  B,  p,  T) 

for  the  location  of  the  solid  carbon  boundary  surface. 


VI.  THE  THERMODYNAMIC  PROPERTIES  OF  GAS  MIXTURES 

The  thermodynamic  properties  of  a  gas  mixture  of  known  composition  are  easily  determined  by 
means  of  well-known  relations.  These  relations  are  discussed  in  detail  in  standard  treatises  on 
thermodynamics.  12/  In  the  present  section,  we  list  without  derivation  the  relations  that  are  con¬ 
veniently  employed  in  the  evaluation  of  the  thermodynamic  properties  of  combustion  gases.  Vie 
include  in  general  form  the  relations  for  taking  into  account  the  effect  of  gas  imp  erfection  on 
these  calculations. 


The  State  Variables 


The  gross  composition  of  the  system  has  been  quantitatively  described  for  the  purpose  of  the 
calculation  of  the  equilibrium  composition  by  the  parameters  qj,  j  =  1,  2,  .  .  .  c,  expressing  the 
composition  of  the  hypothetical  mixture  consisting  of  components  only.  The  normalization  of  the 
qj,  expressed  by  equation  (3.14),  implies  that  qj  is  the  number  of  gram  moles  of  the  j-th  compo¬ 
nent  in  the  hypothetical  mixture  of  components  only  whose  gaseous  phase  is  of  weight  M  ,  where 


M0  -  q;  Mj 

j 


(3-15) 


17/  See  for  example,  MacDougall,  F.  H.,  Thermodynamics  and  Chemistry:  John  Wiley  &  Sons, 
Inc.,  New  York,  3d  ed.,  1939. 

See  also  Beattie,  J.  A.,  Work  cited  in  footnote  13. 
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and  where  Mj  is  the  molecular  weight  of  the  j-th  component.  Explicit  expressions  for  M0  have 
been  given  for  each  of  the  particular  computation  programs.  The  total  number  of  moles  of  gas 
contained  in  a  gas  phase  of  weight  MQ  has  been  shown  to  be 

n  =G_1,  (6.1) 


where  the  general  definition  of  G  is 


s 

G  =  i  +  ( v  i-i)  x£. 

i=c+ 1 


(4.11) 


Explicit  expressions  for  G  have  been  given  for  each  of  the  particular  computation  programs.  We 
have  also  shown  that 


s 

n  =  1  +  >>  ( 1-  V  £)  nj 

i=c  +  l 


(4.16) 


for  a  computation  conducted  in  terms  of  die  number  of  moles  n^  instead  of  the  mole  fractions  xj, 
of  the  constituents  k  =  1,  2,  .  .  .  s.  The  mean  molecular  weight  M  of  the  gas  phase  is  evidently 
given  by 

M  =M0/n.  (6.2} 

We  have  assumed  that  the  gas  phase  is  described  by  a  law  in  one  or  the  other  of  the  alterna¬ 
tive  forms 

p/P  -  -rf—  <P(p,T),  (6.3) 

Mo 

,  n  RT 

p/P  = -rj-  ¥(<PT).  (6.4) 

For  an  ideal  gas,  <p  =  1  and  \y  =  1.  The  first  expression  determines  the  density  as  a  function  of 
pressure,  temperature,  and  composition.  The  second  expression  determines  the  pressure  as  a 
function  of  density,  temperature,  and  composition.  The  molar  volume  v  and  the  density  are  re¬ 
lated  by 

v=M/p=M0/nP.  (6.5) 

The  Thermodynamic  Properties  of  Ideal  Gas  Mixtures 


The  enthalpy  and  energy  of  a  .system  are  undefined  to  the  extent  of  an  arbitrary  additive  con¬ 
stant.  This  fact  is  of  no  consequence  in  practical  applications,  as  one  always  wishes  to  know 
the  difference  of  enthalpy  (or  energy)  between  two  different  states,  and  this  difference  is  unam¬ 
biguously  defined.  The  enthalpy  or  energy  of  an  equilibrium  mixture  of  combustion  gases  must 
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therefore  be  evaluated  relative  to  some  convenient  state  taken  as  standard.  Evidently,  in  cal¬ 
culating  the  difference  of  enthalpy  or  energy  between  two  states,  the  functions  for  each  state 
must  be  evaluated  relative  to  the  same  standard  state. 

As  standard  state  for  the  equilibrium  mixture  of  combustion  gases,  we  shall  usually  employ 
the  stoichiometrically  equivalent  mixture  of  elements,  each  element  being  in  its  customary  stand¬ 
ard  state.  The  standard  state  of  carbon  is  taken  to  be  solid  graphite  at  0°  K.,  and  that  for  hydro¬ 
gen,  oxygen,  or  nitrogen  to  be  hypothetical  ideal  gas  at  O0  K. 

We  define  die  relative  molar  enthalpy  H  of  a  gas  mixture  as  the  difference  between  the  en¬ 
thalpy  of  1  mole  of  the  mixture  and  die  enthalpy  of  the  stoichiometrically  equivalent  mixture  of 
elements  in  dieir  standard  states.  For  an  ideal  gas  mixture,  we  write 

s  s 

H*  (p,T)  =  xk  (H°  -  Eg)k  +  ^  xk  (AfEg)k,  (6.6) 

k=l  k=l 


where  the  asterisk  is  employed  to  specify  that  the  quantity  is  evaluated  for  an  ideal  gas  mixture. 
In  this  expression,  (H°  -  Eg)k  is  the  change  in  enthalpy  when  1  mole  of  the  k-th  constituent  is 
taken  from  its  standard  state  to  a  state  of  temperature  T  and  sufficiently  low  pressure  such  that 
the  ideal  gas  law  and  the  Gibbs-Dalton  law  are  applicable,  and  (AfEg)k  is  the  change  of  energy 
(equal  to  the  enthalpy  at  0°  K.)  accompanying  the  formation  of  1  mole  of  the  k-th  constituent  in 
its  standard  state  from  the  elements  in  their  standard  states.  The  quantity  (H°  -  Eg)k  is  a  func¬ 
tion  of  the  temperature  only  that  can  be  calculated  by  statistical  methods  from  the  spectroscopi¬ 
cally  observed  energy  levels  of  die  k-th  constituent.  (AfEg)k  is  a  constant  for  given  k  that  can 
be  deduced  from  appropriate  colorimetric  experiments.  The  relative  molar  enthalpy  H*  of  the  gas 
mixture  is  a  function  of  both  temperature  and  pressure,  as  the  mole  fractions  xk  are  functions  of 
both  temperature  and  pressure. 

For  an  ideal  gas  mixture,  the  relative  molar  energy  E  of  the  mixture  is  given  by 

s  s 

E*  (P,T)  =  ^  xk  <E°  *  Eg)k  +  ^  xk  <AfEo)k>  (6'7> 

k=t  k  - 1 

where  (E°  -  Eg)k  is  the  change  in  energy  when  1  mole  of  die  k-th  constituent  is  taken  from  its 
standard  state  to  a  state  of  temperature  T  and  sufficiently  low  pressure  so  that  the  ideal  gas 
laws  are  applicable.  For  ideal  gas  mixtures,  we  also  have  the  well-known  relation 

H*  =E*  +  pv  =E*  +  RT,  (6.8) 

where  v  is  the  molar  volume  and  R  is  the  gas  constant. 

The  molar  heat  capacity  at  constant  pressure  Cp*  and  die  molar  heat  capacity  at  constant 
volume  Cv*  can  be  computed  by  the  relations 

s 

Cp*  =^>  xk  Cpg, 
k  =  l 

(6.9) 

s 

Xk  Cvg, 

k=l 
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where  Cpg  and  Cvg  are  the  heat  capacities  at  constant  pressure  and  at  constant  volume,  respec¬ 
tively,  of  the  k-th  constituent  at  temperature  T  and  pressure  sufficiently  low  so  that  the  ideal  gas 
laws  are  applicable.  For  an  ideal  gas  mixture,  we  also  have  the  well-known  relation 

Cp*  =Cv*  +  R.  (6.10) 

The  adiabatic  exponent  Y*  can  be  calculated  from  the  relations 


Y*  =Cp*  +  _  Cp* 

Cv*  Cv*  Cp*  -  li 


(6.11) 


Cpg  and  Cvg  are  functions  of  temperature  only  that  can  be  calculated  by  statistical  methods  from 
the  spectroscopically  observed  energy  levels  of  the  k-th  constituent. 


The  molar  entropy  of  an  ideal  gas  mixture  is  given  by 


S*  (p,T)  -  xk  Sg  -  R  xk  log  xk  -  R  log  p, 
k=l  k=l 


(6.12) 


where  Sg  is  the  entropy  of  the  k-th  constituent  in  the  hypothetical  ideal  gas  at  unit  pressure  and 
temperature  T  and  can  be  calculated  by  statistical  methods  from  the  spectroscopically  determined 
energy  levels  of  the  k-th  constituent. 

In  the  consideration  of  combustion  processes,  it  is  frequently  desirable  to  refer  the  extensive 
thermodynamic  properties  to  a  fixed  weight  of  the  system.  The  specific  enthalpy  h*  is  the  en¬ 
thalpy  per  unit  weight  and  is  given  by 

1  n 

h*  (p,T)  =-  H*  (p,T)  =  ~  H*  (p,T),  (6.13) 

the  specific  energy  e*  by 

1  n 

e*  (p,T)  =  -  E*  (p,T)  =  —  E*  (p,T),  (6.14) 

M  M0 

and  the  specific  entropy  s*  by 

**  <P.D  -  i  S*  (p,T)  -  i.  S>  (p,T).  (6.15) 

The  thermodynamic  properties  of  the  gas  mixture,  assumed  to  be  ideal,  may  be  considered  to 
be  functions  of  the  temperature  and  density.  The  dependence  of  the  properties  considered  on  the 
pressure  (or  density)  is  explicit  only  in  the  case  of  the  entropy.  In  the  case  of  all  of  the  other 
properties,  the  dependence,  on  pressure  (or  density)  is  implicit,  arising  from  the  dependence  of  the 
equilibrium  compositions  on  the  pressure  (or  density).  The  relations  given  can  be  employed  with¬ 
out  modification  for  computation  of  the  enthalpy,  energy  content,  and  heat  capacities  of  the  equi¬ 
librium  mixture  at  specified  density  and  temperature.  The  entropy  is  given  by  the  relation 

s  s 

RT  P 

S*  (  P  ,T)  =  xk  Sg  -  R  xk  log  xj<  -  R  log  “jj  (6. 16) 

k=l  k=l 
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The  Thermodynamic  Properties  of  Imperfect  Gas  Mixtures 


Beattieiiv  has  comprehensively  discussed  the  thermodynamics  of  imperfect  gas  mixtures  that 
contains  a  systematic  derivation  of  the  general  relations  required  to  take  account  of  the  effect  of 
gas  imperfection  on  the  calculation  of  the  thermodynamic  properties.  In  the  present  section  we 
assemble  those  relations  of  interest  in  the  consideration  of  combustion  gases. 

If  die  equation  of  state  is  of  the  form 

Pv  =RT  cp(p,T),  (6.17) 

where  v  =M/  p  is  the  molar  volume,  one  obtains  the  following  relations 

o  rP  dcp  dp 

H(p,  T)  =  H*(p,T)  -  RT 2  /  -J-  JL, 

o  i  p 

S(p , T)  =  S*(p,T)  -  R  /P  [<P  -  1  +  T  ]  —  ,  (6-18) 

o  o  i  p 

*  rP  r  02CP  Q qp  dp 

Cp(p,T>  3  Cp(p,T).  RT 

for  the  relative  molar  enthalpy  H,  molar  entropy  S,  and  molar  heat  capacity  Cp  of  the  real  gas  mix¬ 
ture  at  temperature  T  and  pressure  p.  The  quantities  designated  by  an  asterisk  are  evaluated  as 
described  in  the  preceding  section.  The  integrals  can  be  evaluated  explicitly  for  particular  forms 
of  the  equation  of  state.  The  relative  molar  energy  content  is  evaluated  from  the  relation 

E(p,T)  =H(p,T)  -  RTcp,  (6.19) 


and  the  heat  capacity  at  constant  volume,  if  required,  may  be  evaluated  with  the  aid  of  the  well- 
known  relation 


S-cv  =  T  Wr)v  [Wr 


(6.20) 


If  the  equation  of  state  is  of  the  form 


E  =  H,(P>X)> 


one  obtains  the  following  relations, 


E(p,T)  =  E*  ( P,T)  -  RT2  / 


P  0t|/  dp 

o  df  ~  ’ 


(6.21) 


rP  r  3^  1  dP 

S(  P,T)  =  S*  ( p,  T)  -  R  /  0-1+  T  ]  r— 

O  O  I  p 


(6.22) 


*  ,P  r  S2  V  T  dp 

C„(P,T)-Ct(p,T)-RT  /o  [T— -2 


Ji/  Beattie,  J.  A. ,  Work  cited  in  footnote  13. 
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for  the  relative  molar  energy  content  E,  molar  entropy  S,  and  molar  heat  capacity  Cv  at  tempera¬ 
ture  T  and  density  P  .  The  quantities  designated  by  an  asterisk  are  evaluated  as  described  in 
die  preceding  section.  The  integrals  can  be  evaluated  explicitly  for  particular  forms  of  the  equa¬ 
tion  of  state.  The  relative  molar  entropy  can  be  evaluated  from  the  relation 

H(P,T)  =  E(P,T)  +  RT\]/,  (6.23) 

and  the  heat  capacity  at  constant  pressure,  if  required,  can  be  evaluated  with  the  aid  of  equation 

(6.20). 

VII.  DESCRIPTION  OF  A  GENERALIZED  TABLE  OF  THE  EQUILIBRIUM  COMPOSITION 
AND  THERMODYNAMIC  PROPERTIESOF  COMBUSTION  GASES 


The  Explosives  and  Physical  Sciences  Division  of  the  Bureau  of  Mines  is  calculating  a  mas¬ 
sive  table  that  will  list  the  equilibrium  composition  and  thermodynamic  properties  at  equilibrium, 
over  an  extended  range  of  temperature  and  pressure,  of  gas  mixtures  containing  compounds  of  car¬ 
bon,  hydrogen,  nitrogen,  and  oxygen.  The  methods  of  computation  described  in  this  report  have 
been  employed  in  the  tabulation.  12/ 

The  table  includes  all  possible  combinations  of  the  four  elements,  the  spacing  of  calculated 
points  being  small  enough  to  permit  accurate  interpolation  in  the  table.  The  temperature  ranges 
from  1,000°  to  5,000°  K.  and  the  pressure  from  0.1  atmosphere  to  pressures  great  enough  to  in¬ 
clude  those  of  interior  ballistics.  The  table  is  the  result  of  nearly  a  million  individual 
computations. 

The  results  of  the  individual  computations  are  recorded  on  punched  cards.  Because  of  the 
size  of  the  project,  publication  in  tabular  form  of  the  complete  table  is  impractical.  However,  the 
file  of  punched  cards  constitutes  a  reference  table  with  all  of  the  computational  results  in  easily 
accessible  form.  This  reference  table  can  readily  be  employed  in  the  consideration  of  specific 
problems.  The  results  of  such  applications  will  be  published  in  a  series  of  reports  to  be  issued 
by  this  laboratory. 

It  is  anticipated  that  the  table  will  have  a  wide  range  of  application  to  the  special  interests 
of  many  organizations.  The  cooperation  of  interested  organizations  is  invited  so  that  the  greatest 
possible  utilization  of  the  table  may  be  achieved.  The  Bureau  will  welcome  suggestions  for  ap¬ 
plying  the  table  to  specific  problems. 

The  state  is  uniquely  determined  by  specification  of  two  state  variables,  for  whichpurpose 
we  select  the  temperature  and  pressure.  In  the  table,  the  temperature  takes  on  the  values 

T(°K.)  =  1,000  (100)  1,600  (200)  5,000, 

where  the  interval  spacing  of  the  argument  of  a  table  is  given  in  parentheses,  and  the  lower  and 
upper  limits  for  which  the  interval  applies  are  given  to  the  left  and  right,  respectively,  of  the 

19/  A  preliminary  announcement  of  this  table  appeared  in  Chem.  Eng.  News,  vol.  27,  Sept.  5, 
1949,  p.  2540-2541. 
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number  in  parentheses.  This  statement  is  to  be  read  "The  temperature,  in  degrees  Kelvin,  ranges 
from  1,000  to  1,600  with  an  interval  of  100  and  from  1,600  to  5,000  with  an  interval  of  200.’’ 

A  closer  spacing  is  necessary  at  low  than  at  high  temperatures  to  permit  accurate  interpola¬ 
tion  with  respect  to  temperature  in  the  table.  It  has  been  found  by  means  of  sample  calculations 
that  interpolation  with  respect  to  the  pressure  is  facilitated  by  the  use  of  a  logarithmic  scale. 
In  the  table,  the  pressure  ranges  from  0.1  to  100  atmospheres  in  the  following  way: 

log  p  =-1.0  (0.2)  2.0. ■ 

It  is  intended  subsequently  to  extend  the  portion  of  the  table  that  is  of  interest  in  considering  gun 
propellants  to  pressures  high  enough  to  permit  its  application  to  the  calculations  of  interior  bal¬ 
listics.  In  this  portion  of  the  table  the  temperature  and  density  will  be  taken  as  the  independent 
state  variables. 

In  selecting  the  table  arguments  that  define  the  gross  composition,  we  have  employed,  for  the 
heterogeneous  region  containing  solid  carbon,  the  quantities  D  and  E  given  by  equations  (3.3).  In 
the  table,  these  variables  take  on  the  values 

D  =0(0.1)  0.9 
E  =0  (0.1)  1.0. 

For  the  homogeneous  part  of  the  table  we  employ  the  parameters  A,  B,  and  C,  defined  by  equa¬ 
tions  (3.2).  In  the  table,  the  variables  A,  B,  and  C  take  on  the  following  values: 

A  =0(0.1)  0.9 

B  =0  (0.1)  0.5,  1.0 

C  =0. 1  (0.05)  0.9  subject  to  C 

For  each  composition  defined  as  indicated  and  for  each  temperature  and  pressure,  the  compo¬ 
sition  at  thermodynamic  equilibrium  is  calculated  by  the  numerical  methods  that  have  been  de¬ 
scribed  in  this  report.  We  have  assumed  that  the  mixture  at  equilibrium  may  contain  the  following 
gaseous  constituents:  CO,  CO2,  H2,  H2O,  N2,  C>2,  OH,  NO,  H,  N,  O,  NH3,  CH4.  This  selection 
neglects  a  number  of  additional  possible  constituents,  and  it  was  dictated  by  the  storage  capacity 
of  the  computer.  However,  it  is  believed  that  the  constituents  of  major  importance  have  been  con¬ 
sidered.  In  exceptional  cases,  where  it  is  desired  to  consider  the  existence  of  additional  constit¬ 
uents,  the  present  work  will  provide  the  starting  basis  for  a  more  detailed  calculation,  which 
should  be  very  simple. 

The  mole  fraction  of  each  constituent  of  the  gas  at  equilibrium  is  recorded  in  punched-card 
form.  In  addition,  for  the  heterogeneous  mixtures  the  punched  cards  record  the  value  of  the  de¬ 
pendent  gross  composition  parameter,  which  defines  the  location  of  the  carbon  boundary  surface. 


The  results  of  the  calculation  of  the  equilibrium  composition  are  employed  in  the  usual  man¬ 
ner  to  calculate  the  density,  energy  content,  heat  content  (enthalpy),  and  entropy  of  the  equilib¬ 
rium  mixture.  The  table  of  thermodynamic  properties  may  be  considered  to  be  a  Mollier  chart,  in 
punched-card  form,  for  each  of  the  different  mixtures  of  the  four  elements  considered. 

In  constructing  the  table  we  employed  the  values  of  the  thermodynamic  properties  of  the  indi¬ 
vidual  constituents  recommended  by  the  National  Bureau  of  Standards..??/  For  pressures  below 
100  atmospheres  we  employed  the  ideal  gas  equation  of  state. 

Only  brief  reference  can  be  made  here  to  some  of  the  applications  anticipated  for  the  table. 
The  flame  temperatures  of  a  variety  of  fuels  burning  under  a  wide  range  of  conditions  can  readily 
be  estimated.??/  Hottel  and  Eberhardt??/  have  discussed  the  application  of  Mollier  diagrams  to 
performance  studies  of  internal  combustion  engines,  and  Sucrow??/  summarizes  their  application 
to  performance  studies  of  jet-propulsion  and  gas-turbine  engines.  The  table  will  be  useful  in  pre¬ 
dicting  the  operating  conditions  for  various  gas-synthesis  processes  that  will  produce  a  product 
gas  of  desired  characteristics. 


20/  Tables  of  Selected  Values  of  Chemical  Thermodynamic  Properties:  Issued  at  irregular  inter¬ 
vals  by  the  National  Bureau  of  Standards,  Washington,  D.  C. 

21/  Lewis,  B.,  and  von  Elbe,  G. ,  Combustion,  Flames,  and  Explosion  in  Gases:  Cambridge  Uni¬ 
versity  Press,  Cambridge,  1938,  Chap.  XIX. 

22/  Hottel,  H.  C,,  and  Eberhardt,  J.  E.,  Mollier  Diagrams  for  the  Internal-Combustion  Engine: 
Chem.  Rev.,  vol.  21,  Dec.  19  37,  pp.  439-460. 

23/  Sucrow,  M.  J.,  Principles  of  Jet  Propulsion  and  Gas  Turbines:  John  Wiley  &  Sons,  Inc.,  New 
York,  1947,  Chap.  III. 
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